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ABSTRACT 

This  study  presents  the  results  of  an  experiment  in  objective 
analysis  of  oceanographic  data  for  a  limited  area.   The  objective 
analysis  is  designed  to  provide  a  reliable  operational  system  for 
tactical  use  in  coastal  waters.   It  is  shown  that  this  approach  makes 
it  possible  to  obtain  a  very  detailed  analysis  with  good  vertical 
consistency  and  that  only  a  relatively  small  amount  of  highly  accurate 
data  is  required.   The  procedure  includes  a  polynomial  interpolation 
and  a  dynamical  interpolation.   The  dynamical  interpolation  which  is 
introduced  involves  the  analytical  solution  of  a  finite  differential 
equation  by  means  of  a  rotational  computing  molecule.   The  procedure 
requires  only  a  small  computer  and  little  computer  time.   This  method 
will  provide  a  basis  for  short-time  forecasts  of  oceanographic  para- 
meters using  only  small  computer  centers  or  even  time  sharing  systems. 


TABLE  OF  CONTENTS 

I .  INTRODUCTION  7 

A.  BACKGROUND  9 

B.  PURPOSE  AND  SCOPE  OF  THIS  STUDY  11 

II .  APPROACH  15 

A.  POLYNOMIAL  INTERPOLATION  15 

B.  DYNAMICAL  INTERPOLATION  17 

1 .  Analytical  Mode  1  18 

2 .  The  Numerical  Model  27 

III .  DESCRIPTION  OF  DATA  34 

IV.  DISCUSSION  OF  RESULTS  36 

A.  COMPARISON  OF  OBJECTIVE  AND  SUBJECTIVE  ANALYSIS  36 

B.  INTERPRETATION  OF  RESULTS  IN  TERMS  OF  THE 

ANALYTICAL  MODEL 41 

C.  COMPARISON  OF  OBJECTIVE  ANALYSIS  RESULTS  FOR 

DIFFERENT  LEVELS  43 

D.  OBJECTIVE  ANALYSIS  OF  SURFACT  TEMPERATURE,  SALINITY 

AND  SIGMA-T  FOR  LOW  LATITUDES 46 

V.  CONCLUSIONS  AND  SUGGESTION  FOR  FUTURE  RESEARCH  54 

COMPUTER  OUTPUT  56 

COMPUTER  PROGRAM 73 

BIBLIOGRAPHY  81 

INITIAL  DISTRIBUTION  LIST 83 

FORM  DD  1473  85 


LIST  OF  FIGURES 


Figure  Page 

1  Example  of  grid  used  by  the  Joint  Numerical  Weather 
Prediction  Unit  (after  Cressman  1959) 8 

2  California  coast  Cruise  Track  showing  station  locations 
(after  Shepard  1970) 12 

3  Cruise  Track  showing  station  locations  along  Region  II 

of  Caribbean  coast  of  Colombia  (after  Sanchez  1969)....      13 

4  Sequence  of  rotation  of  the  computing  molecule  around 

a  grid  point  with  no  value  U(K,L) , 32 

5  A  complete  rotation  of  the  computing  molecule 33 

6  Subjective  analysis  of  sea  surface  temperature  for  28 
April-4  May  1970  (after  Shepard  1970) 37 

7  Objective  analysis  of  sea  surface  temperature  for  April- 
May  1970,  showing  wind  field.   Arrowy  are  the 

vec  tor 38 

8  Comparison  of  subjective  and  objective  analyses  of  sea 
surface  temperature  for  California  coast.   28  April- 

4  May  1970 39 

9  Objective   analysis   of   temperature,    California  coast,    at 

the    10-meter   level    for   28  April-5  May   1970 45 

10  Superposition  of  main   features    from  surface   and   10-meter 
level   analyses,    California  coast.      Thin   lines    correspond 

to   10-meter   level 47 

11  Objective    analysis   of  sea  surface    temperature,    California 
coast,    for  September   1969 48 

12  Objective   analysis   of    temperature    at   the    10-meter   level, 
California  coast,    for   September   1969 »  49 

13  Geographic  position  of   the   grid   for   the   Low  Latitude 
Analysis 50 

14  Objective   analysis   of   sea  surface    temperature    for  Low 

Lati  tude 51 

15  Objective   analysis   of  sea  surface   salinity   for    low 

lati  tude 52 

16  Objf  ctive   analysis   of   sea  surface   sigma-T   for   low 

lat:  tude 53 


AKNOWLEDGEMENT 

It  is  time  to  thank  the  United  States  Navy  for  making  it  possible 
for  the  author  to  be  at  the  Naval  Postgraduate  School,  the  School's 
faculty  for  giving  him  knowledge  and  means  to  carry  out  this  work 
especially  to  his  Co  advisor  Assistant  Professor  Kenneth  L.  Davidson, 
Department  of  Meteorology,  whose  patience,  understanding  and  collab- 
oration made  possible  the  solution  of  so  many  theoretical,  semantical 
as  well  as  personal  problems  that  the  author  had  to  face  during  the 
development  of  this  work. 

The  author  wants  to  thank  Professor  Dale  F.  Leipper,  Chairman  of 
the  Oceanography  department  and  Co  advisor  of  this  work,  for  his  in- 
valuable advice  and  continuous  cooperation  not  only  during  the  com- 
pletion stage  of  this  work  but  also  during  all  the  stages  of  its 
development  including  its  planning;  also,  Mr,  Stevens  P.  Tucker  who 
kindly  supplied  all  the  data  for  the  California  coast  used  in  this 
work;  the  W.  R.  Church  Computer  Center  staff  because  without  their 
continuous  help  and  highly  efficient  performance  this  work  wouldn't 
be  possible;  the  Education  Media  Department  staff  who  made  possible  a 
neat  presentation  of  the  drawings . 

The  author  also  wants  to  express  his  personal  recognition  to 
Admiral  Jaime  Parra  Ramirez,  Colombian  Navy,  to  whom  the  Oceanography 
in  Colombia  owes  its  existence  and  whose  interest  in  this  field  may 
have  influenced  the  general  approach  and  objectives  of  this  work. 

Finally,  it  is  in  order  to  state  the  author's  deepest  recognition 
to  Lieutenant  Colonel  Fred  B.  Tomlinson,  U.  S.  Army,  from  the  Letterman 


General  Hospital  surgery  department  who  performed  an  extremely  success- 
ful operation  on   the   author's    three-month  old  daughter  during  critical 
stages   of   this  work  making  possible   to   finish  it  with  clear  mind  and 
without  personal  worries. 


I.   INTRODUCTION 

Analysis  of  oceanographic  data  has  been  traditionally  a  very  slow 
and  essentially  a  subjective  process  due  to  the  time  involved  in  an 
oceanographic  expedition  and  the  small  amount  of  data  obtained.   In- 
creasing the  number  of  systematic  observations  in  time  and  space 
requires  faster  analysis  procedures  and  objective  methods.   Such  methods 
will  provide  information  for  regions  of  sparce  data.   The  analysis  pro- 
cedure should  also  standarize  sequential  analyses  so  that  the  results 
may  be  compared. 

Objective  analysis  techniques  are  currently  used  on  an  operational 
basis  by  the  Navy  and  other  government  agencies  concerned  with  atmos- 
pheric phenomena.   Background  information  available  for  formulating 
an  objective  analysis  scheme  for  oceanographic  data  is  primarily  that 
from  inves-tigations  and  from  experience  with  schemes  for  meteorological 
data.   However,  there  are  significant  differences  between  both  the 
type  and  distribution  of  meteorological  and  oceanographic  data  which 
have  to  be  considered. 

The  primary  purpose  of  most  objective  analysis  schemes  is  to 
transform  data  from  observations  at  irregularly  spaced  points  into  data 
at  points  within  a  regularly  arranged  grid  (grid  points).   A  subjective 
method  for  doing  this  would  be  to  have  an  analyst  manually  (subjectively) 
draw  isolines  and  then  interpolate  (objectively/subjectively)  to  obtain 
values  at  the  appropriate  grid  points.   An  example  of  an  operational 
grid  is  given  in  Figure  1. 
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FIGURE  1.   Example  of  the  grid  used  by  the  Joint  Numberical  Weather 
Prediction  Unit  (After  Cressman  1959) 


Procedures  used  in  objective  analysis  have  to  be  systematic  and 
quantitative  so  that  they  can  be  performed  by  numerical  computations 
on  a  computer.   Computer  processing  is  by  its  very  nature,  objective. 
In  practice,  once  the  observations  are  relayed  from  reporting  stations 
all  procedures  leading  to  an  analyzed  chart  should  be  completed  within 
a  minimum  of  manual  intervention. 

A.   BACKGROUND 

Research  and  the  subsequent  development  of  objective  analysis 
methods  took  place  in  response  to  the  need  to  prepare  large  numbers 
of  meteorological  observations  representing  the  atmosphere,  both  in 
the  horizontal  and  vertical,  for  numerical  weather  prediction  procedures 
Numerical  weather  prediction  requires,  in  addition  to  values  at  reg- 
ularly spaced  grid  points,  editing  (removal  of  erroneous  observations) 
of  these  data  and  checking  for  consistency  and  continuity  in  time  and 
space.   These  requirements  would  also  apply  to  objective  analysis 
schemes  for  oceanographic  data. 

Once  a  grid  of  equally  spaced  values,  checked  for  possible  errors, 
is  generated  by  an  appropriate  objective  analysis  scheme  further  work 
can  be  performed  readily  on  a  computer.   These  steps  also  could  be 
performed  in  subjective  analysis  but  the  time  required  to  obtain  the 
completeness  achieved  by  an  objective  analysis  scheme  would  be  out  of 
proportion  to  the  time  needed  for  the  computer  analysis. 

An  early  publication  on  the  subject  of  objective  analysis  of  geo- 
physical data  was  by  Panofsky  (1949).   Various  approaches  for  per- 
forming objective  analysis  with  meteorological  data  were  presented 
during  successive  years  by  Gilchrist  (1954),  Bergthors  r>n   and  Doos 
(1955),  Jomson  (1957),  Cressman  (1959),  and  Smith  (1962).   The  latter 


references  describe  versions  of  operational  objective  analysis  schemes 
presently  used  in  the  principal  weather  centers  of  the  world.   However, 
most  of  the  currently  used  objective  analysis  schemes  are  based  on  the 
"Operational  objective  analysis  system"  described  by  Cressman  (1959). 

The  general  approach  in  these  schemes  is  a  very  straightforward 
corrective  approach  in  which  an  initial  (preliminary)  field  is  modified 
by  available  observations  and  the  result  is  then  used  as  preliminary 
field  which  is  again  modified  with  the  observations  but  with  more 
resolution.   The  repeating  corrections  are  called  scans.   This  process 
may  be  repeated  many  times  but  in  general  satisfactory  results  are 
obtained  with  the  first  or  second  scan. 

In  the  corrective  approach,  an  analytical  model  defines  the  cri- 
terian  for  checking  the  validity  of  observed  values  and  determining 
values  at  grid  points.   Determination  of  grid  point  values  from  obser- 
vations at  nearby1  stations  is  achieved  by  some  form  of  interpolation. 
One  method_of  interpolating  for  a  scalar  field  is  to  represent  the 
contours  of  isolines  surrounding  the  grid  point  as  a  function  of 
position.   This  may  be  accomplished  by  means  of  a  surface  least  squares 
polynomial  interpolation.   The  order  or  degree  of  the  polynomial  in  the 
least  square  interpolation  depends  on  how  detailed  the  final  field 

should  be; 

"The  degree  of  smoothing  will  depend  upon  the  nature  of  the  fore- 
cast and  the  particular  model  selected Smoothing  to  eliminate 

unwanted,  though  real,  wave  lengths  has  been  taken  to  fall  within 
the  province  of  the  forecast  rather  than  that  of  the  analysis. 
(Johnson  1957) 


1  Nearby  is  normally  considered  to  be  within  one  grid  distance. 
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Oceanographic  analyses  have  been  made  only  as  a  related  factor  in 
weather  analysis.   Climato logical  forecasts  of  ocean  surface  temperature 
anomalies  have  been  made  by  Adem  (1969),  Jacob  (1967),  and  Adem  (1970) 
using  the  equation  for  the  conservation  of  thermal  energy. 

■   The  type  of  data  expected  in  an  oceanographic  survey  introduces 
the  problem  of  holes  between  transects  and  therefore  the  need  of  some 
type  of  dynamical  interpolation.   Typical  cruise  tracks  appear  in 
Figures  2  and  3.   Provisions  for  lack  of  data  for  meteorological 
analysis  of  a  given  region  when  good  data  is  available  at  the  boundary 
have  been  studied  by  Thompson  (1961),  Richardson  (1961),  Smith  (1961), 
and  Smith  (1962).   The  method  generally  used  is  to  "advect  good  data 
from  outside  the  region  (or  "hole")"  (Smith  1962). 

B.   PURPOSE  AND  SCOPE  OF  THIS  STUDY 

The  purpose  of  this  study  is  to  develop  an  objective  analysis 
scheme  for  oceanographic  data  which  will  meet  operational  requirements 
described  in  the  following  paragraphs. 

The  high  cost  or  operating  oceanographic  ships  at  sea  is  a  primary 
reason  to  develop  methods  to  complete  instant  detailed  analysis.   Such 
an  analysis  would  show  interesting  features  based  on  early  observations 
but  not  considered  during  the  initial  planning  which  would  permit 
scheduling  of  complementary  observations  during  the  same  operation  at 
sea.   Furthermore,  operation  requirements  of  submarine  and  antisubmarine 
warfare  have  shown  the  need  of  faster  and  more  accurate  short  term 
forecasts  of  oceanographic  parameters  which  could  be  produced  through 
the  use  of  numerical  procedures. 

The  scope  is  to  develop  a  scheme  which  would  resolve  fine  details. 
Small  thermal  inhomo gene i ties  have  great  effect  in  sound  propagation 
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especially  when  relatively  high  frequencies  are  used.   Those  fine 
details  cannot  be  observed  nor  forecasted  with  large  scale  schemes  and 
therefore,  this  study  considers  small  scale  objective  analysis  methods 
for  oceanographic  data.   The  method  will  include  some  type  of  dynamical 
interpolation.  The  scheme  will  be  evaluated  with  respect  to  its  use 
for  tactical  purposes  and  as  a  "feed  back"  for  planning  during  ocean- 
ographic  surveys. 
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II.   APPROACH 


Data  obtained  from  the  standard  ooeanographic  stations  differ  from 
the  meteorological  data  in  the  nnmber  of  stations,  the  accuraey  of 
observations  and  in  the  fact  that  the  former  nearly  always  has  to  be 
considered  synoptic1  for  long  time  steps  over  small  areas.   In  general, 
an  initial  field  will  not  be  available  in  ooeanographic  analysis  nor 
will  there  be  enough  data  to  obtain  a  complete  interpolated  field. 
Therefore,  new  procedures  have  to  be  formulated  in  the  development  of 
the  objective  analysis  for  ooeanographic  data.  Techniques  described 
above  for  meteorological  analysis  will  be  applicable  only  in  a  relatively 

small  part  in  this  particular  problem. 

._  i ioTM>d  Is  =   numerical  scheme  in  which  a 

The  objective  anarysrs  wv^"^  tB  -  "J 

least  squares  surface  polynomial  is  used  for  interpolating  and  smoothing 
the  available  data  to  nearby  grid  points.   Then  a  dynamical  inter- 
polation is  used  to  compute  values  at  those  grid  points  which  are  in 
"holes"  or  regions  of  no  data.   The  dynamical  interpolation  is  achieved 
by  means  of  a  rotational  molecule. 

A.   POLYNOMIAL  INTERPOLATION 

The  first  step  in  the  scheme  is  to  reduce  the  original  data  to 
neighboring2  grid  points  using  surface  least  squares  polynomial 


1  Synoptic  data  represent  observations  taken  simultaneously. 

2  Neighboring  refers  to  data  within  one  grid  distance. 
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interpolations.  A  second  degree  polynomial  is  used  when  siX  or  .ore 

daCa  points  are  available  in  a  given  square  around  the  grid  point,  and 

first  degree  polynomial  when  only  £our  or  five  points  are  available. 

A  weighted  gradient  propagation  is  used  when  only  two  or  three  data 

points  are  available.  When  only  one  observation  is  available  in  the 

entire  square  around  the  grid  point  some  criteria  about  the  distance 

U   established  and  if  this  distance  is,  say,  less  than  five  miles  the 

•11  hP  to  eive  the  same  value  to  the  grid  point, 
best  approximation  will  be  to  give 

otherwise  this  value  is  neglected. 

^e  polynomial  solved  for  the  second  degree  polynomial  interpolation 


is: 


a„  ♦  4.x«  a.  xy  *  am  x'+  <K,f  a.,  f    « 


The 


first  degree  interpolation  formula  is: 

~am  ia„x+  a„xy  +  a„y 


(2) 


and   the  weig 


ighted  gradient  propagation  formula  is: 


w*\**g&g!?&   " 


Uhere  ]rf=J.O        *»«•  «*«  "^   °f  ^  ""  P°  *  C°nSldSred' 
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The  matrix  obtained  in  the  first  two  cases  is  solved  for   doo 
since  the  grid  point  is  at  location  0,0    ,  using  the  subroutine  SIMQ 
from  the  computer  library  of  the  Naval  Postgraduate  School. 

This  scheme  for  obtaining  values  for  some  grid  points  enables  the 

maximum  use  of  data. 

In  meteorological  analysis  the  regression  has  been  used  and  also 
the  correction  method.   In  the  latter,  it  is  assumed  that  the  difference 
between  the  preliminary  value  and  the  observed  value  at  a  given  station 
is  equal  to  the  difference  between  the  derived  value  and  the  preliminary 
value  at  the  grid  point.   Another  approach  in  the  correction  method  is 
to  assume  that  the  gradient  of  the  preliminary  field  at  the  grid  point 
is  representative  for  all  the  area  between  grid  points  and  the  station 
(gradient  propagation).   Then  the  correction  is  a  function  of  the 

,  -    .        i  _  f  n  ■--'  ->—    T«  f-Vio««a  ni*of*pdu'JL6S  a 
distance  between  grid  points  anu  station,   in  cnese  p.o_e_ui.„ 

statistically  determined  weight  function  (weighted  gradient  propagation) 

is  required . 

In  all  interpolation  methods,  the  distance  between  observations  and 

a  grid  point  becomes  a  crucial  factor;  for  example, 

"Quite  often  it  is  not  possible  to  get  a  reasonable  analysis  only 
by  means  of  interpolation  between  synoptic  observations.   It  is 
quite  clear  that  the  distance  between  observations  must  be  small 
compared  with  the  size  of  the  system  to  be  analyzed." 
(Bergthorsson  and  Doos  1955) 

Therefore,  distance  limitations  introduce  the  need  for  a  dynamical 
interpolation  whenever  "holes"  are  present. 

B.   DYNAMICAL  INTERPOLATION 

The  initial  step  in  developing  a  suitable  dynamic  interpolation  is 
simplification  of  the  basic  conservation  equations  for  a  scalar  quantity 
in  a  fluid  medium.   In  order  to  apply  a  simplified  equations  to  a  set  of 
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regular  horizontally  and  vertically  distributed  data,  certain  assump- 
tions are  required  about  the  flow  characteristics.   These  assumptions 
have  to  be  based  on  conditions  observed  or  generally  accepted  for  an 
oceanographic  environment.   The  simplification  of  these  equations,  the 
assumptions  involved,  and  the  numerical  procedures  used  to  apply  the 
final  form  of  the  equation  are  discussed  in  the  following  paragraphs. 

1.  Analytical  Model 

Consider  the  analysis  of  a  conservative  scalar  field.   The 
available  data  will  represent  quasi-synoptic  observations  provided  by 
standard  stations  in  a  relatively  small  area  of  the  ocean.   Therefore, 
some  of  the  following  information  could  be  available  for  each  station: 

(1)  Temperature 

(2)  Salinity 

(3)  Wind  velocity 

(4)  Dynamic  height  anomaly 
(1)  Sigma-t 

(6)  Sound  velocity 
at  observation  depths  and  at  standard  depths. 

An  analytical  model  for  analysis  of  temperature,  salinity  and 
sound  velocity  would  be  based  on  the  following  steady  state  conservation 
equation  for  a  non  divergent  fluid: 


^■x     1U         "?*     2* 


K 


111 


(*) 
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where 


u,  v,  and  w  are  the  components  of  the  velocity  vector 
K  is  the  eddy  viscosity  coefficient 

Studies  have  shown  that  the  following  approximation  is  possible: 


2. 


K  ^]  *  \a/% 


^2 


(5) 


where 


"The  condition  for  neglecting  the  horizontal  terms  is  that  they 
approximately  cancel  each  other.   Since  those  terms  have  been 
shown  to  be  of  the  same  order  of  magnitude  and  opposite  signs, 
this  condition  is  possible."  (Overstreet  and  Rattray  1969) 

Assuming  condition  (5)  the  conservation  equation  becomes: 


9x     u^% 


-*\   r 


"WO 


j4 


~\  /  \  i 


(6) 


Assuming  isotropy  and  horizontal  homogeneity  the  two  dimensional 
Equation  (6)  becomes: 


Z  /  \  "\ 


(7) 


or 


2u+  vdn]  =  2b+!& 


(8) 


ier    Uto 
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where     r  ~      -—- 

u 

P  can  be  computed  for  each  grid  point  if  the  direction  of  the  current 
at  each  grid  point  is  known.   However,  unless  direct  measurements  of 
the  current  have  been  made  near  each  grid  point  this  is  not  generally 
the  case.   Therefore,  the  next  best  solution  is  to  compute  the  current 
for  each  station  using  available  data  and  extrapolate  or  interpolate 
it  to  the  nearest  grid  points  assuming  that  the  current  is  uniform  in 

the  area. 

The  current  consists  of  four  components 

V-  w  +  y+vL+Vr  do) 


where 

Yj      is    the   drift   current  due    to  wind  action 

U      is-the   geostrophic   current  computed   from  the  density   field 
6 

Vi       is    the  barotropic  current 

K      is    the   tide  current 
If   a  strong  barotropic   current,     VL      ,    is   present   in   the   region   to  be 
analyzed,    climatological  data  would  have    to  be   considered   for   its 

definition. 

Tidal  currents  may  be   significant  because  oceanographic   data 
is   only  quasi-synoptic   as   noted  before.      However,    currents   considered 
are   generally   time   averaged   for   the  entire   time  of   all   the  observations. 
Therefore,    observations  would   represent   a  time   lapse   of   two  or   three 
days   or  at   least  24  hours.      It  is   possible,    therefore,    to  neglect   tidal 
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treats   -  *.  average.     However,  when  the  ohservations   «  «* 

synoptic    tidal     currents  must  be  considered,   i.   .. 

„Tldal  currents   are  simply  --""STSSl'SSS  are'peri 
forces.      Since   tidal  forces  vh ^         lraportant  in  shallow 
odical  and  these  currents   are  gen «       *         d  as  part  o£   the  ad- 
water,   tidal  currents  -U^o tjbe  e^^  Forecasting,  Asvffirs 
vection  problem.        <>uceQ™ 
manual  series,   v.    5,    p.    *■).) 

Hereafter  we  neglect  tidal  and  harotropic  currents  so  the 

horlaontal  components  of   the  current  vector  become-. 

a  -  Uj  +  ^ 


(id 


nr'TTd*^ 


p.  3&JM 

similar  formulae  have  heen  suggested  hy  Arthur  (1966,  and  Adem 

(1970).       -  . , 

+»  velocity  components   can  be   coated   -  the  *.«  avail- 
able  in   a  standard  oceanographic   station. 

u      .  ,f,   current     we   examine   possible   relation 
Considering   the  drift   current,   v 

♦-v,0   cnrface  wind  velocity 
hetween  it  and   the  obsereved  wind.     Denoting  the  surfac 

as  va,    the   following  formula  apply   for  Va  along  the  Y-axrs: 


&  =  ve 


(12) 
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where 


T)    -  JMLS 


and 


V-  JKOJgjL±Vd 

T)  =  i*JJ<i 


fcr   \V.\  >6/rn/M<. 


(J)        is    the    latitude 

y  .... 

,i   ,»  ,.,.-..■-  -•  ~  «<■><-  alone  the  wind  airectK 

For  the  general  case,  when  the  i-ctxj.o  j-o  not  a±.ong  uw  

Equations  (12)  becomes  (Adem  1970), 


(13) 


d    ^^  ^  J 

Where  $     is  the  angle  between  the  surface  current  and  the  surface 
wind.      #    is  positive  when  the  current  is  to  the  right  of  the 

surface  wind. 

The  geostrophic  current  (   ^  ,  dj    >,  between  two  stations 
may  be  computed  from  the  dynamic  height  anomaly  using  the  Ekman 
approximation: 
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(14) 

where     f     is  the  Coriolis  parameter 

AT)  is  the  geopotential  anomaly 

Assuming  that  the  density  increases  steadily  with  depth,  the  relative 
geostrophic  current  between  two  stations  A  and  B,  can  be  expressed  as 

ar-%-^\-A^  (15) 


where  L  is  the  distance  between  A  and  B 

A  further  simplification  of  Equation  (11)  may  be  achieved 
assuming  that  the  geostrophic  current  has  no  contribution  to  the 
advection  of  density   anamolies; 

"The  major  contribution  to  advection  (of  heat)  comes  from  the 
wind  drift  current  since  the  slop  current  is  assumed  to  be 
geostrophic,  or  flowing  rarallel  to  the  isotherms.   That  is, 
there  is  no  advection  (of  heat)  by  water  flowing  parallel  to  the 
isotherms;  only  cross  section  flow  causes  a  change  xn  the  tem- 
perture  at  a  point,"  (Ocean  Thermal  Structure  Forecast,  SP-105, 
ASWEPS  manual  series,  v.  5,  p.  33) 

If  this  approximation  can  be  accepted,  Equation  (11)  becomes, 


(16) 
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7>=   JU 

Therefore,  many  degrees  of  approximation  may  be  used  to  compute 
P  and  they  can  be  summarized  in  three  cases: 

a.  currents  at  the  station  are  known.   If  the  true  current  has 
been  measured  at  each  station,  P  may  be  computed  directly  or  if  the 
current  has  been  computed,  P  may  be  computed  from  Equation  (11).   The 
same  P  may  be  used  to  nearby  grid  points  provided  the  changes  in  current 
direction  are  not  too  large  relative  to  the  grid  scale. 

b.  The  wind  velocity  is  kn™™  but  the  true  current  is  not 
k.own  for  each  station.   In  this  case  only  the  drift  current  can  be 
considered  and  Equation  (16)  must  be  used  to  obtain  the  direction  of  the 
current  for  the  grid  points  using  the  wind  direction  nearby.   Also, 

in  this  case  empirical  consideration  may  be  used  for  allowing  ] 
drift  current.   These  empirical  considerations  are:   When  slope  and 
gradient  currents  (due  to  mass  transport  by  the  wind)  have  been  set  up, 
they  affect  the  direction  of  the  original  wind  driven  current,  and  when 
a  steady  state  has  been  reached,  the  deflection  of  the  surface  from 
the  down  wind  direction  will  be  reduced  to  near  20  degrees  at  low 
wind  speeds  or  to  less  than  12  degrees  at  high  wind  speeds  (Summary  of 
Technical  Reports  of  the  National  Research  Committee,  v.  6A) ,  (Ocean 
Thermal  Structure  Forecasting,  SP-105,  ASWEPS  manual  series,  v.  5), 

(Hubert  1964). 

c.   Only  the  wind  field  is  known  but  ot-he.r  factor  may  be 
important.   This  situation  may  arise  when  the  coast  or  bottom  topography 
would  invalidate  the  assumption  that  the  current  is  at  a  constant  angle 
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to  the  wind  direction.   For  this  case  it  is  possible  to  use  the  follow 
ing  relation  to  obtain  as  estimate  of  the  direction  of  the  current  at 
each  grid  point. 


K^  V 


(Kitaigorodskii  1961) 


and  therefore 


J  Vi      -   COS  TAN  T    =    RATIO 

K 

where  V  is  the  fluctuating  velocity  vector. 

However,  sufficient  surrounding  data  are  required  to  apply  a  numerical 
solution  for  obtaining  these  variables.   In  order  to  compute    ]V± 
a  P  value,  from  the  wind  field  is  assumed  for  each  grid  point  as  in 
b.  or  for  the  whole  grid  if  the  wind  field  is  more  or  less  uniform. 
Using  the  assumed  P,  values  of  U/K  are  computed  for  all  possible  grid 
points  using  Equation  (8).   The  numerical  sheme  to  perform  this  will 
be  explained  in  the  next  section. 

Defining  y  U    7   +    W  4 

T~-  tA     l<i 

and  assuming  that  eddies  present  in  the  area  may  introduce  a  great 
range  of  directions  in  the  fluctuating  velocity  and  given  that   IX- 
is  a  constant,  the  maximum  value  of  U/K  should  correspond  to  a  minimum 
value  of  W/K       •   Also,  if   ^   is  small  enough,  it  is  possible  the 


approximation 


IV]    -v   U    maximum      —     RATIO 
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This  approximation  can  be  made  if  it  is  assumed  that  the  mean  current 
is  thermohaline  and  therefore  flows  parallel  to  isotherms  so  it  does 
not  contribution  in  the  advection  of  heat.i.  e. 

\/SVl  +  V'  VT  =    V'T 
if    V  1  VT 

Therefore,    only   the   fluctuating   component   is   obtained  when       V 
and         VT  are   computed   from  a  field  of   temperature    to  solve   for 

V  .      For  each  grid  point   at  which  U/K  is    computed   aOJ/K  can  be 

obtained  using   the   formula 


ijk     r.  n2      ,a2t^ 


¥-\m-d)W-x* 


(17) 


K 


and   the  new  P   =     —rf  -      ,. 


Using  this  U/K  and  the  new  P  Equation  (8)  can  be  solved  for  T  at  next 
grid  point  if  enough  nearby  grid  values  exist  to  use  the  rotational 
computing  molecule  presented  in  the  numerical  model  (seccion  2). 

The  procedure  above  can  be  repeated  many  times  but  it  has  been 
observed  that  usually  the  first  approximation  is  satisfactory.   Because 
the  quadratic  equation  (17)  has  two  solutions,  the  wind  field  is  used 
in  the  selection  of  the  correct  solution.   This  choice  can  be  objective 
and  included  in  the  computer  program. 
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2.   The  Numerical  Model 

Knowing  P,  U/K  can  be  computed  whenever  the  parameter  is 
Known  in  four  neighboring  grid  points  and  then,  using  this  value  of 
U/K  and  the  corresponding  P  the  next  unknown  grid  point  parameter  can 
be  computed.   The  problem  reduces  to  solving  the  differential  equations 


u 


2Q+vlfJ 


(8) 


ui2cj  =  HJ 


(5A) 


The  corresponding  Finite 

the  same  scheme . 

The  following  is  a  description 


Differential  Equations  may  be  solved  using 

of  a  direct  solution  of   the 


par 


tial  differential  equations  by  means  of  a  rotating  compu 


tational 


mo 


lecule  in  such  a  way 


that  a  solution  is  possible  even  though  data  are 


scarce. 


The  horizontal  continuity  equation,  Equation  (8),  may  be 
approximated  by  the  following  computational  molecules: 
(1)  a. 
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(2)    a 


,    2(K 


~K*XR      l-  J 


b.- 


z 


(4)   a.    J  n  -  (2^flz(^UH)-2.0*Z(^l)fZ(^U'Oj 

^2.0^2(K-1JL)~2CK-2JL)-  M^JiK-l^-liK-^L) 

k*xscl 

b.  7(xi)  =  ?tK0-  -4£ — [z(x-<^)-z(x-2^-z)] 


(5)    a.     Zfc0=    ^^^^-^^^^-^"2>^^?; 


K*XR    l  J 
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b.     7(K,L)--    J  (KtL)  -  (U\lT*f  \z(K,L-Z-)-Z(K-2,L)\ 


"    Z(K*>)  =  %(*>*)-£ 


u 


0*K+XSC 


2(K*U-i)-2(x-tjL-i)i 


(7)   a.  2p(KtL)^  40*Z(Ki1,L-l)-Z(Ki2>L)-Z(KJLr2)-'2(Ki^Lr2) 


_  (tt*nrt 


*^  r  ?  A*  K*Xh 


2.0*K*XR 


2(**W-2(Xjt-'Q 


(8)    a. 


or  "I 

K*xscL  J 


(9)    a 


2.o*K*yscl  J 
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I„  each  case   the  computing  molecule  has  been  divided  in  two  parts   in 
order   to   give  more  flexibility  to   the  computation  and  be  able  to  compute 
the  Lapacean  and  the  gradient  separately.      In  order  to  be  able  to  use 
all  molecules,    the  X-axis  scale,  XSC,  must  be  equal   to   the  Y-axis 
scale,  YSC.     If  the  preceeding  is  not   the  case,    local  distortion  may  be 
introduced.     Also  the  above  procedures   are  not  iterative  and,    therefore, 
stability  problem  will  not   arise. 
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FIRGURE  4.   Sequence  of  rotation  of  the  computing  mole  lie  around 
a  grid  point  with  no  value  U(K,L). 
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FIGURE 


5,      A  complete   rotation  of    the   computing  molecule 
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III.   DESCRIPTION  OF  DATA 

In  order  to  verify  the  proposed  scheme,  data  were  selected  from 
areas  with  good  data  coverage  and  well  known  physical  features. 

One  area  considered  was  the  California  coast  between  San  Francisco 
Bay  and  Monterey  Bay  and  the  corresponding  data  were  obtained  by  the 
AGOR1  cruises  of  September  1969  and  28  April  -  4  May  1970.   These  data 
correspond  to  standard  oceanographic  stations  (Figure  2).   However, 
not  all  the  information  was  readily  available  due  to  the  fact  that  the 
data  were  reduced  for  specific  uses.   The  data  available  were; 

a.  Wind 

b.  depth 

c .  temperature 

d.  Sound  velocity 

e.  Beam  transmittance 

at  observation  depths  only.  Under  these  circumstances  the  general 
procedure,  Equation  (11),  doesn't  apply  and  therefore,  the  modified 
procedure,  c,  was  used. 

A  second  area  was  considered  in  order  to  verify  the  objective 
analysis  procedures  at  low  latitudes.   This  second  set  of  data  was 
obtained  by  the  Colombian  Navy's  Division  of  oceanography  in  the 
Caribbean  sea  off  the  Colombian  coast.   In  this  case  all  the  data  from 


1  Cruises  by  ships  of  the  U.  S.  Oceanographic  Office  under  inves- 
tigation plans  by  the  Naval  Postgraduate  School. 
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a  standard  oceanographic  station  were  available.  However,  the  same 
simplified  prooedure,  c,  was  used  in  order  to  obtain  results  com- 
parable  to  those  for  the  California  coasts  case. 
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IV.   DISCUSSION  OF  RESULTS 

A.   COMPARISON  OF  OBJECTIVE  AND  SUBJECTIVE  ANALYSIS 

A  subjective  analysis,  Shepard  (1970),  of  sea  surface  temperature 
for  the  area  along  the  California  coast  appears  in  Figure  6.  Shepard 
performed  the  analysis  in  order  to  make  a  comparison  of  oceanographic 
parameters  during  upwelling.   The  objective  analysis  of  the  same  data 

appears  in  Figure  7. 

For  the  objective  analysis,  Figure  7,  the  grid  was  rotated  47  degrees 
in  order  to  follow  the  general  trend  of  the  coast.   Because  the  mean 
wind  was  generally  along  the  X-axis  during  the  cruise,  an  initial 
value  of  P  equal  -0.27  or  a  current  20  degrees  to  the  right  of  the 
X-axis  was  selected. 

A  comparison  clearly  shows  greater  complexity  and  more  details  in 
the  objective  analysis  than  the  subjective  analysis.   The  obvious 
question  that  arises  is;  do  these  details  represent  artificial  additions 
or  real  features  with  physical  meaning?  Also,  do  they  represent  features 
expected  with  the  inherent  assumption  of  the  objective  scheme? 

It  may  be  stated  that  subjective  analysis  is  a  linear  interpolation 
which  normally  depends  only  on  the  available  data  for  one  oceanographic 
variable  at  a  time.   This  generally  excludes  consideration  of  currents. 
Representations  in  intermediate  regions  with  no  data  are  generally 
smoothed  contours.   The  resulting  scalar  field  (as  depicted  in  Figure  6) 
is  very  familiar  in  the  oceanographic  work  and  a  departure  from  this 
picture  is  often  viewed  with  suspicion.   Adverse  appraisal  of  a  more 
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•„   „^t-    t-nt-allv   iustified  as   stated  in  the 
detailed  representation   is  not   totally   jusciiLeu 

following  passage: 

STSSSSi  structur^becomes    apparent  but  no, |  -^  *«£»- 

by  the  data  that   the  analysis  becomes  difficult  (Ocean  Tnerm 

Structure  Forecasting,  ASUEPS  manual  series,  v.    5,   p.   "J. 

Other  features  present  in   the  objective  and  absent  in  the  subjective 

analysis   are   the  circular  patterns  or  pools  of  relative  warmer  or  cooler 

water.     This   type  of  temperature  distribution  is   intuitively  logical 

as  most  swimmers  can  tell  and  physically  may  be  expected.      Subjective 

analysis  procedures  are  not  suitable   for  detecting  them  because   they 

do  not  normally  involve  considerations  of   temperature  along  with 

advection  currents   and  diffusion  terms. 

_   ,,       -,    -i.__.--j       TTcnnllv   a   single    line 
Mnjj,-,.  v,ove  seldom  been   tuny   delineated,      usuanj    *  ft  B    .        -,    ^ 

2Z£^  -r-r^rrf  t:izi  iu/=£5_sr£r 

r_b^:__c_."Wa«  of"  Technical  Reports  of  the  National 
Defense  Research  Committee,  v.  6A,  p.  7 J) 
The  importance  of  including  the  influence  of  the  wind  action,  i.  e. 
wind  currents,  within  the  scalar  analysis  is  emphasized  further  in  the 
following  physical  arguments. 

Wind  and  wave  action  have  a  net  result  of  mass  transport  of  water 
and  therefore  both  are  capable  of  producing  zones  of  convergence  or 
divergence  accompanied  by  a  rising  or  lowering  of  the  thermocline  or 
may  be  responsible  for  water  of  one  temperature  being  carried  over 
water  of  another.   (Sugary  of  Technical  Reprots  of  the  National 
Defense  Research  Conmittee,  v.  6A) .  For  example, 

"The  effect  of  the  wind  in  maintaining  or  deepening  an  isothermal 
surla   layer  is  complicated  by  the  fact  that  the  <**"«£  * 
current  near  the  surface  which  may  remove  the  mixed  ater  from 
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r      -    „*    *r   i.    formed."      (Summary  of   Technical 
KS--3  r^f  Defend  tSS  CoJ.ee,  v.   6A,   p.  *) 

B.      INTERPRETATION  OF   RESULTS   IN  TERMS  OF  THE  ANALYTICAL  MODEL. 

Additional  interpretation  o£   the   features  within  the  objective 
analysis   are  possible  by  considering   the  basio  processes   included  in 
the  analytical  oodel.     The  physical  description  of   the  analytical  model, 
Equation   (4),  used  in  the  objective  analysis  scheme  is   that  when  a  steady 
state  is  reached  the  advection  of   the  scalar  property  is  equal   to  the 
eddy  diffusion  of   the  scalar  in  the  medium.     With  this  physical  inter- 
pretation, we  can  explain  certain  features  of   the  circular  distribution 
o£  properties  when  these  features   are   frozen  in  time   (stationary). 

These   features    are: 

1.  Circular  patterns, 

2.  Tongue  type  patterns, 

3.  Gradients 

4.  Siz.e  of  closed  isolines. 

With  respect  to  the  circular  patterns,  heating  and  cooling  of  the 
sea  surface  is  not  uniform  for  the  entire  layer  of  the  analysis  and 
therefore  in  any  event  a  local  circulation  will  be  generated  in  order 
to  obtain  the  equilibrium  of  mass.  Warm  water  will  flow  from  the 
source  toward  the  periphery  in  all  directions  and  the  Coriolis  accel- 
eration will  deflect  the  paths  to  the  right  (in  the  Northern  Hemisphere) 
until  they  become  circular.   This  situation  will  exist  as  long  as  a 
disequilibrium  in  mass  exists  and  will  become  a  steady  state  when 
advection  out  of  the  circle  is  equal  to  the  inflow  of  warm  water,  a 
similar  situation  is  reached  in  a  cold  pool  where  warm  water  from  the 
periphery  is  going  to  the  center  of  the  pool  followin  a  circular  path. 
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Decaying  drift  currents  will  also  produce  the  same  effect  because  once 
the  wind  ceases  the  drift  current  continues  to  move  for  a  certain  time 
under  its  inertia  while  being  acted  upon  by  the  Coriolis  acceleration 
and  while  being  damped  by  the  diffusion  term. 

With  respect  to  tongue  type  features,  the  oceanographic  literature 
has  abundant  examples  of  tongue-like  diffusion  of  properties  such  as 
salinity  and  temperature.  For  example,  it  has  been  observed  that  in 
areas  of  convergence,  convection  cells,  the  surface  water  moves  down- 
wind at  a  considerably  more  rapid  rate  than  in  areas  of  divergence  and 
therefore  the  flow  of  water  in  response  to  the  wind  should  be  uneven. 

With  respect  to  the  gradients,  because  vertical  mass  transport 
occurs  in  the  convective  layer,  it  is  clear  that  the  assumption  of  a 
horizontally  homogeneous  mixed  layer  above  the  thermocline  is  an  over- 

l     •    „l,       -.„_  AS  ^rro      maw      T-.A      DOSSible      ill 

simplified  concept  and  relatively  high  gradients  ..a,  __  pc 
short  horizontal  distances.   This  is  especially  true  in  the  area 
examined  which  is  in  coastal  waters,  where  an  irregular  wind  field 
may  have  existed  and  where  upwelling  is  known  to  exist.   The  analytical 
model  was,  of  course,  a  non  divergent  and  this  may  not  have  been  the 
case.   Therefore,  greater  than  normal  gradients  inside  the  inertial 
circles  may  reflect  over  compensation  for  the  upwellinS  or  downwelling 
effect  in  the  temperature.  Although,  this  is,  perhaps,  an  unrealistic 
result,  it  could  be  smoothed  out  by  computing  the  amount  of  this  effect 
using  the  grid  points  over  nearby  coastal  regions  (U.  S.  Fleet  Numerical 
Weather  Facility.   Technical  memo  No.  5,  January  1965).   On  the  other 
hand,  artificial  intensification  of  the  gradient  may  still  be  used  to 
identify  convergence  and  divergence  zones  by  allowing  a  maximum  gradient 
in  a  given  region  depending  of  local  conditions. 
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With  respect  to  the  size  of  closed  isolines,  the  radius  of  the 
inertlal  circle  for  the  type  of  circulation  mentioned  above  varies 
with  the  velocity  and  the  sine  of  the  latitude:   For  a  current  of 
average  velocity  located  in  mid-latitudes  this  radius  is  about  15  miles 
which  can  be  greater  if  the  circulation  Is  influenced  by  a  horizontal 
component  of  gravitational  force  or  other  lateral  force  (Summary  of 
Technical  Reprots  of  the  National  Defense  Research  Committee,  v.  6A) . 
It  is  interesting  to  note  that  15  miles  is  the  order  of  magnitude 
observed  in  the  circular  features  of  the  objective  analysis  under 

consideration. 

Therefore,  it  seems  possible,  with  this  analytical  model,  to  relate 
circular  and  elliptical  features  to  inertlal  currents  due  to  differ- 
ential heating  and  cooling,  convergences,  divergences,  decaying  wind 
drift  currents;  tongue  like  features  to  normal  diffusion  processes; 
and  intensified  gradients  inside  circular  patterns  to  strong  upwelling 
or  downwelHng.  However,  gradients  are  represented  very  well  within 
the  analytical  model.   Reasons  for  the  latter  feature  are  however,  known. 

C.   COMPARISON  OF  OBJECTIVE  ANALYSIS  RESULTS  FOR  DIFFERENT  LEVELS 

In  order  to  verify  vertical  consistency,  a  mixed  layer  of  at  least 
10  meters  was  assumed  and  an  objective  analysis  of  the  10-meter  level 
was  made  and  is  shown  in  Figure  9.  However,  the  wind  field  was  not 
uniform  and,  there  was  not  enough  information  about  fetch  and  duration 
with  available  wind  observations.  Therefore,  we  do  not  expect  to  obtain 
an  exact  correspondence  of  the  features  of  the  two  layers  over  the 
entire  grid.  Nevertheless,  vertical  continuity  is  observed  when  the 
main  features  of  the  two  analyses  are  superimposed  in  ^igure  10. 
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Maximum  and  minimum  temperature  regions  are  marked  with  letters  A 
through  K  in  Figure  7  and  apparent  corresponding  features  in  Figure  9 
have  been  marked  with  same  letters  with  a  prime  sign.   The  same  regions 
are  shown  in  Figure  10  but  delineated  by  only  one  or  two  isotherms  in 
order  to  simplify  the  resulting  picture. 

The  comparison  of  thermal  structure  at  the  two  levels  shows  an 
almost  identical  correspondence  between  regions 

A  and  A' 
B  and  B' 
D  and  D1 
E  and  E1 
F  and  F1 
H  and  H1 
I  and  I1 
J  and  J1 
This  may  be  expected  from  good  vertical  mixing  and  more  or  less  constant 
advection  in  the  entire  layer. 

Regions        C  and  C1 
K  and  K* 
G  and  Gf 
have  horizontal  relative  displacement  from  5  to  12  miles.   This  may  be 
due  to  a  different  velocity  field  for  the  two  levels.   Past  history 
(i.  e.  two  or  three  weeks  earlier)  would  show  how  these  advective 
distortions  were  generated.   Similar  results  may  be  obtained  for  the 
September  data  (Figures  11  and  12). 

Considering  that  each  analysis  is  completely  independent  from  the 
other,  this  comparison  shows  convincing  verification  «   vertical 

consistence. 

44 


u 

o 

m 

r-( 

CI) 

> 

<u 

1-t 

u 

<u 

•u 

Q) 

B 

i 

o 

r-t 

<u 

js 

•U 

4J 

Ctf 

»* 

■U 

CO 

« 

o 

o 

CO 

•H 

g 

o 

M-l 

•rl 

i-l 

CtJ 

o 

*\ 

CD 

^ 

3 

i-> 

ctJ 

^ 

<u 

a 

e 

CD 

■U 

• 

<4-l 

o 

o 

r^ 

o^ 

co 

r-) 

•H 

W 

>, 

>^ 

CCj 

r-< 

s 

cd 

C 

LPl 

cO 

a) 

> 

l-l 

•H 

•H 

4J 

h 

CJ 

CX 

CD 

< 

•i-l 

*o 

CO 

O 

CN 

CT> 

B 

a 

H 


45 


D.   OBJECTIVE  ANALYSES  OF  SURFACE  TEMPERATURE,  SALINITY  AND  SIGMA-T 

FOR  LOW  LATITUDES 

Objective  analyses  of  temperature,  salinity  and  sigma-t  for  a  region 
with  mean  latitude  of  10°  north  appear  in  Figures  14,  15,  and  16.   The 
main  reason  for  considering  these  analyses  for  a  region  not  very  well 
studied  from  the  oceanographic  point  of  view  is  to  compare  the  effect 
of  latitude  in  the  size  of  those  features  interpreted  to  be  inertial 
circles  in  the  preceding  discussions. 

The  region  of  California  coast  considered  has  a  mean  latitude  of 
37°  north  and  feature  A  in  Figure  7  appears  to  be  representative  one  for 
this  region.   It  is  observed  that  its  mean  radius  is  about  8  nautical 
miles.   The  radius  of  inertial  circles  is  inversely  proportional  to  the 
sine  of  the  latitude  and  therefore  it  would  be  expected  that  the  radius 
of  the  inertial  circle  be  increased  by  a  factor  or  approximately  3.4 
when  translated  from  latitude  37°N.  to  luV  a  corresponding  feature 
in  Figure  14  was  marked  A"  has  mean  radius  of  about  35  nautical  miles. 
Therefore,  the  test  for  latitude  effect  seems  satisfactory  from  this 
gross  comparison.   However,  more  precise  computations  would  be  required 
before  a  conclusive  result  may  be  obtained. 

The  region  off  the  Colombian  coast,  considered  in  the  low  latitude 
study  (Figure  13),  receives  the  run-off  of  three  main  rivers:   the 
Atrato  in  the  south-west  region  of  the  grid  and  the  Sinu  and  Magdalena 
rivers  in  the  east.   Features  L,M  and  N  in  Figure  15  correspond  to  places 
where  the  run  off  occurs.   The  patterns  observed  at  these  points  coin- 
cide very  well  with  what  may  be  expected  and  are  very  well  correlated 
with  the  sigma-t  distribution,  L' ,  M'  and  N'  in  Figure  16. 
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V.   CONCLUSIONS  AND  SUGGESTION  FOR  FUTURE  RESEARCH 

A.   CONCLUSIONS 

A  small  scale,  high  resolution  and  low  computing  time  objective 
analysis  model  has  been  developed.  More  tests  and  field  verification 
are  required  to  show  all  of  its  advantages  and  shortcomings  but  tests 
described  and  interpreted  in  this  study  suggest  that  it  has  great 
possibilities  to  become  a  basic  step  in  forecasting  oceanographic 
parameters  on  a  tactical  unit  level.   One  desirable  feature  of  the 
results  was  the  detail  obtained  in  contrast  with  that  obtained  in  a 
subjective  analysis.   The  required  input  and  simplicity  enables  it  to 
be  used  to  obtain  an  almost  instantaneous  analysis  during  an  oceanographic 
survey  and  therefore  to  investigate  interesting  features  not  foreseen 
in  planning  stages.   Therefore,  it  could  be  applied  for  the  purposes  of 
reducing  costs  and  increasing  scientific  results  of  an  oceanographic 

expedition. 

The  procedures  described  need  to  be  refined  in  order  to  allow  for 
convergence  or  divergence  whenever  these  phenomena  become  important. 
However,  by  understanding  how  the  analytical  model  (section  B)  handles 
these  features,  meaningful  interpretations  can  be  made  even  for  these 
features  in  the  results.   The  data  should  be  corrected  for  diurnal 
effects  but  there  is  not  a  satisfactory  procedure  for  this  (Technical 
report,  ASWEPS  No.  13).   An  alternative  solution,  observations  at  the 
same  hour  of  the  day,  may  be  undertaken  by  XBT  drops.   This  is  perhaps, 
unfeasible  at  most  oceanographic  stations  due  to  operational  costs. 
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B.      SUGGESTIONS   FOR  FUTURE   RESEARCH 

A  suggestion  for  future   investigation  would  be   to  integrate   this 
scheme  into   a  forecast  procedure.     The   following  steps  could  be  applied 
using  the  present  form  of   the  scheme: 

a.  Objective  analysis  using  horizontal  grid  at  standard  depths. 

b.  Conversion  to  a  three-dimensional  grid  using  the  vertical  grid 

as    a  complementary   tool. 

c.  Repetition  of  procedures   a.    and  b.    at  equally  spaced  time  steps. 

d.  Determination  of  Lagrangian  time  rate  of  change   for  identificable. 
structures    (assuming  that  some  of   them  are  persistent  in  relation  to   the 
time  step  selected  as   seems   to  be   the  case  in  the  neighborhood  of 

constant  currents). 

e.     Determination  of  changes   in    fte  heat  budget  for  each  character- 

,     ..     .       ____=!..  „~^,,>-£'i   nsrsmptcrs   such   as 
istic  structure  and  Relating)   it   to  easily  measured  pa 

amount  of  clouds   in  each  standard  level,   sun's  declination   (season)   and 

hour  of   the-  day. 

£.  Measurement  of  advection  by  means  of  changes  in. the  structures 

in  point  d.  and  (relating)  it  to  permanent  currents,  wind  speed,  direction 

and  duration. 

g.  Forecast  for  the  next  time  step. 

h.  Forecast  verification  and  obtention  of  corrective  cefficients. 
in  the  delineation  of  the  steps  leading  to  a  forecast  it  has  been 
assumed  that  some  features  may  have  a  relative  long  duration  or  per- 
sistence but  this  situation  must  be  investigated  for  each  application. 
If  it  is  proved,  identification  of  persistent  or  semipersistent  features, 
their  changes  in  time  and  parameters  affecting  them  may  be  a  simple 
approach  for  detailed  forecasts,  even  for  handling  tr  ,sient  thermocline. 
now  practically  impossible  to  forecast. 
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PROGRAM  OUTPUT 

MESH       MINIMUM  MAXIMUM       HEANLAT  28  April    -   4  May 

n"  a    .  nn  14-000  37.UJU 

EPSILON 


2  6.L0O  14.u00 


ROTATION    PARA,    GRID       INT       INITIAL       INCREMENT       EPSILUN 

POSITICN    CF    THE    GRID    REFERENCE    POINT 
LATITUDE  LONGITUDE 

37  42.200      124  6. ODD 

OCEANOGRAPHIC  STATIONS 


LATITUDE 


LCNGITUDE         DEPTH  TEMP   SOUND  VELOC 

36  40.50   121  53.40  0.0  IJ.JJ  itllllS  otc 

\%  »:§8  til  gig  8:1  ifcli  i«5:X8  8:8 

3?  49.80       121  54.50  0.0  10."  1492.10  0.0 

36  52. 9u       1^1  I4*  15  q    Q  1L07  1491.20  0.0 

36  54.10       121  54.40  ^.~  Jj;g2  1491.10  0.0 

|g  55.70      121  34. 4u  O.g  J    ^  1488.70  0.0 

36  53.  5u      121  iMk  0   0  9*91  1487.40  0.0 

ll  51*00       121  59.10  g.O  go  l490.30                 6.50 

36  56.20       Ui  |9. 30  g.O  iQ*53  1489.40  12.90 

36  54.80       121  59. K-  0.0  lt.||  1487.50  27.70 

36  47.10       121  l\m%'±  or  10.03  1487. 8'.)  29.^ 

36  45.30       121  58. Jy  0.0  ig.v  i487.90  22.30 

36  42.60      121  59.40  g.O  ig.u  1489.10  2C.20 

36  41.20      121  59.20  0.0  lO.hQ  1489.70  24.20 

36  39.40      121  59.20  D.U  1490.50  14.40 

3*  38.20       121  ^v.lL  u.u  iY    37  1492.40                 5.6C 

11  36.00       122  1.60  g.O  11.37  JJ-J.  j^.gO 

36  33.60       122  4.60  O.u  11.  1487.60  64.90 

36  31.30      122  ,I»28  n   n  10^39  1488.30  y5l#?^ 

36  35.30      122  10.70  0.0  iu.a*  1488,70  48.50 

36  39.70       122  15.20  0.0  IJ^J  1A87.40  63.00 

36  42.-70       122  18.  7y  0.0  g.J|  1487o7C  71.70 

36  47°.C0       122  14.70  g.O  lij.i|  1A89.40  79.10 

|g  50.00      122  10.20  O.C  l^.||  148B.30  54.0C 

36  51.60      122  9.70  0.0  it«|  1489.00  26.  5u 

36  53.30       122  9.20  0.0  ig",66  1490.60                 9.70 

36  54.30       122  7.50  1.00  l^«cg  1A92o70                 5.10 

^a  55.70       122  6,y^  0.0  TT  1491.90                 5.7u 

36  55.90      122  4.80  0.0  11.24  [Xai.U  21.40 

37  0.40  122  14.30  0.0  *•£*  1486.70  22. 5u 
37  3.8^  122  17.90  g.O  *•»-  ! 486.50  23. 00 
37  4.20       122  19. 4y  C.O  |.  1486.30                 7.60 

36  57.40       122  10.  Jv;  C.u  9*38  I486. 3U              36.00 

37  2.70  122  23.50  0.0  J  '  1466.90  33*2g 
37  lc50  122  21.80  0.0  9.-6  J^/#40  26.60 
37  0.1C  122  23.30  O.v  * • ^  1487.90  19.30 
^k  59„4L  122  24.0'-  i*Xu  i  r,  z.  i  1489.60  v.-J 
?6  56130  122  28.60  0.0  Jg.41  {^UllO  18. 40 
§f,  53.20      122  33.40  0.0  !£•«»  1494.4C              65.90 

36  58.40       122  40.51  g.O  !£•  1495.70               62.  90 

37  2.60       122  4le^"  n    n  12*22  1495. OC              ^-".V 
37  10.00       122  .45.80  o.O  l|.|g  1495.3u              63.50 

37  10.00       122  ^0?V  oo  11.7-s  1493.30               6£°Z> 

37  11.00       122  38. 5o  0.0  fflkl  1492.60                 9.5C 

37  10.90      122  ^5.2J  C.0  ll»gj  1469.9,;                9.60 

37  11.00      122  33.20  0.0  10.|J  1489.20              10.2C 

37  U.00       122  29.10  0.0  "•  1487. 9u               16.*0 


37 


10.40      122  26.70 
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37 

10.90 

122 

2  5.20 

1.00 

9.35 

1486.00 

26.80 

37 

15.80 

122 

27.  50 

V  •  V 

10.26 

1489.40 

11.00 

37 

21.  60 

122 

31.7c; 

CO 

10.54 

1490.20 

9.50 

37 

31*90 

122 

33.41- 

9.92 

1494.30 

19.80 

37 

31.60 

122 

35.20 

0.0 

11.12 

1492.30 

4.10 

37 

32.10 

122 

37.10 

CO 

11.65 

1493.60 

5.00 

37 

31.90 

122 

3  8„70 

t.o 

11.35 

1492. 3u 

5.30 

37 

32.2C 

122 

41.30 

0.0 

11.41 

1492.60 

22.10 

37 

31.90 

122 

43.10 

CO 

11.18 

1492.50 

31.80 

37 

31.60 

122 

46.  9C 

0.0 

9.76 

1486.20 

62.30 

37 

32.00 

122 

49.90 

0.0 

10.49 

1469.60 

63.10 

37 

32.30 

122 

54.50 

CO 

10.98 

1491.20 

61.20 

37 

3 1  o  2  G 

122 

53.00 

0.0 

11.88 

1493.90 

68.20 

37 

32.30 

123 

3.80 

0.0 

12.27 

1495.10 

72.60 

37 

32.40 

123 

9.30 

CO 

12.06 

1494.40 

63.80 

37 

36.80 

123 

15.40 

CO 

11.07 

1491.60 

64.70 

37 

42.20 

123 

7.70 

CO 

10.46 

1489.40 

77.20 

37 

47.70 

122 

59.70 

Co 

10.19 

1489.20 

6.10 

37 

52.60 

122 

52.50 

0.0 

9.94 

14  8  8.30 

12.00 

37 

53.2C 

122 

49.20 

O.u 

9.94 

1486.30 

9.8U 

37 

49.00 

122 

48.80 

CO 

9.58 

1467.00 

17.70 

37 

46.  40 

122 

46.90 

0.0 

8.95 

1484.60 

49.60 

37 

43.50 

122 

44.  50 

CO 

10.71 

14  90.30 

16.00 

37 

40.20 

122 

42.10 

CO 

10.81 

1490.00 

19,10 

37 

37.  50 

122 

40.20 

0.0 

10.80 

1490.30 

23.00 

37 

37.30 

123 

36.  60 

0.0 

10.62 

1489.50 

18.10 
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RELATIVE  POSITION  IN  THE  GRID 


118. C91 
113.213 
112.482 

110.691 
108.642 
1C7.631 
106. A3  1 
105.316 
106.722 
103.396 
104.529 
110.541 
111.693 
113.286 
114.419 
115.735 
116.667 
116.914 
117.035 
117. 409 
112.468 
106.799 
102.698 
101.733 
1C1.991 
101. C94 
100.123 
IOC. 3 18 
99.948 
100.619 
92.153 
87.7C5 
86.595 
96.689 
87.093 
87.262 
87.469 
87.599 
87.360 
87.012 
79.341 
71.966 
67.971 
70.150 
71.217 
73.088 
74.104 
76.338 
78.084 
78.536 
73.7C0 
67.17C 
58.712 
57.951 
56.550 
55.825 
54.189 
53.428 
51.577 
49.6  50 
46.924 
45.822 
41.858 
38.788 


34.677 

38.244 

38.926 

40.378 

42.727 

43o370 

44.461 

40.040 

39.C56 

41.940 

41.102 

36.259 

34.856 

32.605 

31.767 

30.539 

29*779 

26.818 

^3.429 

20.4  59 

21.026 

21.399 

21o4ul 

26.671 

31.345 

32.729 

34.  13 'J 

35.855 

37.511 

38.523 

3 1>  j  044 

36. 261 

35.653 

36.510 

33.992 

32.414 

30.583 

29.697 

24.896 

19.978 

19.378 

17.629 

24.195 

26.531 

29.141 

31. COO 

32.236 

34.631 

35.623 

36.840 

38.339 

40.343 

4t»©376 

45.120 

44.351 

43.280 

41.966 

40.710 

38.286 

36.807 

34.325 

31.530 

28.893 

25.74° 


PARAMETER 

11.640 
1C.280 
10.310 
10.620 
11.310 
11.070 
11.0  20 
10.270 
9.910 
10.820 
10.5  30 
9.930 
10.030 
10.09u 
10.430 
10.600 
10.750 
11.370 
11.040 
10.070 
lu.390 
10.410 
9.980 
10.120 
10.530 
10.2  50 
1U.600 
10.86  0 
11.400 
11.240 
9.620 
9.550 
9.490 
11.420 
9.  380 
9.560 
9.7  JC 
9.S5L 
10.410 
10.480 
12.c40 
12.450 
12.2  2L 
12.300 
11.700 
11.210 
10.810 
10.210 
9.8  50 
9.  350 
10.260 
lb. 540 
9.920 
11.120 
11.650 
11.350 
11.410 
11.180 
9.760 
10.490 
10.980 
11.880 
12.270 
1 2  . 1 6C 
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22.247 
32.494 
32.830 
33.170 
34.529 
37.818 
40.754 
44.183 
47.903 
50.913 
20.322 


25. 

188 

33. 

368 

41. 

792 

49. 

340 

51. 

676 

49. 

C45 

48. 

381 

47. 

,805 

46« 

,956 

46c 

,224 

13, 

,147 

11.070 

10.46C 

10.190 

9.940 

9.940 

9.580 

8.950 

lo.710 

10.810 

10.800 

10.62L 
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INTERPOLATED  VALUES  FCR  THE  GRID 

Y  PARAMETER 

h  eoaoo  8.C0000  o«o 

J'.cooco  i*aps?g  8:8 

P4100000  80O0OGO  u.u 

15  88000  8.00000  0.0 

aoIcoooo  i-ssRss  8:8 

Jslocooo  S-88PR8  uo 

IS*gS§88  l:?oS8§  o:co 

ca   r^rrtfj  8.00000  «•? 

104  00000  80  ceo 00  o.c 

11?   COOOO  8.00000  0.0 

"8:88888  p04  frUoo 

letocgco  l;f*8888o  ioZ  62800 

24.00000  16.0   jyu  xg« 

32.00000  IS*S9.S88  olo 

40.00000  }fcSrrt?n  olo 

48.000C0  i£°9R8$l8  CO 

56.00000  ^•^Pno  olo 

8:88888  |:c||§  \£\m 

13:88888  K:M  if-.Ssooo 

KS:888S8  ||:M  \$™ 

1  Of*     ft  1  I60  W'JVW  i  *  / 

128  .COO  i?*9Srtftr  riT< 

elcoooo  li- 88888  8:8 

16.00000  §MSfcR8  olo 

2^00000  »      ,:  li:07000 

SS'.SoMO  IfcfcoSSo  12.06000 

48100000-  Z^.OGOOC  0.0 

&-8S8S8  itlc-ooo  12:22000 

tI'.COGCG  ItoOOCO  12.30000 

80.0C0C0  .  24.C00OG  xfcSlOOC 

S6.CC0C0  SMv?O0  1C^74?S 

128.CC000  ^•CCGCO  u.c 

ft    GGCGG  ^2.wjGvjO  V«X 

16.00000  io    'r^nr  CO 

24.00000  32«y\92«  9  ^547 

32C0C000  32.U.0v0  i|-27000 

40.00000  li'nrrrS  lll43289 

48.00000  32„0^bC0  olo 

56.0CCCO  If* 89888  if. 14965 

64.00000  2§*9xSS8  11      >?377 

72.10000  l?*k8S88  9:05577 

80.0C000  li«S8888  9.63905 

88.C00CC  l?-9Xr28  9:§2117 

104.00000  IHnOCO  1§:l4487 

118:88888  11:88888  12:34712 

licOOOCO  32.00CCO  0.0 
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8.0CCCC 
16.C0CwO 
24.00000 
32.00000 
40.0u00t 
48.00000 
56.0C000 
64.00000 
72.0CCC0 
80.00000 
88.0CQC0 
96.vCUtO 
104.000^-0 
112.CC000 
120.00000 
128.CO0JC 
6.00000 
16.CC000 
24.CC0C0 
32.C000C 
40.000C0 
48.00000 
56.CC00U 
64.00000 
72.C0000 
80. C  C000 
88.0v00u 
96.CC000 
104.00000 
112.CCCuy 
12C.0C0CC 
128.00000 


40.C0C00 

40. COO  00 

4Q  .  u  w  .^u 

40.00000 

40.00000 

40 .000 00 

40.CC00C 

40.00000 

4  3. C 0000 

40.CG000 

40. COO 00 

40.  COO  00 

40.00000 

40.CGOOO 

4C.C0CC0 

40.00000 

48.C0CC0 

48o0CC0G 

4B.C0C00 

48. (COCO 

48.0u00C 

48.00000 

48.CG0OO 

48.  I CO  CO 

48.00000 

48.00000 

48.00000 

48. CCO 00 

48.0C000 

48.00000 

48.  uoO-jO 

48. COO 00 


0.0 

0.0 

0.0 
10.19000 

0.0 
10.26933 
11.43151 
10.54000 
10.26000 

9. 35C00 

9.55000 

0.0 
10.54430 
10.45130 
0.0 
CO 
0.0 
0.0 
0.0 
10.00981 
8.11982 
10.72071 
8.92048 
10.39016 
0.0 
0.0 
0.0 
0.0 

1L020CC 
12.20851 
0.0 
CO 
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U/K 
0.6C0 
G.6GC 
0.GC6 
0.6C0 
0.6CO 
G.6CC 
0.6CC 
0.591 
C.6CC 
■C.033 
0.600 
0.6CG 
C.6CC 
■0.374 
0.6C0 
0.600 
0.155 
0.192 
0.600 
0.6CC 
O06CU 
0.6C0 
C.6C0 
0.547 
0.546 
Co  600 
0.6C0 
0.600 
•0.064 
0.6CC 
0.6CC 
0.6CC 
0.600 
0  •  6  G  C 
C.6CC 
C.6D0 
C.6C0 
G.6CC 


V/K 
0.0 

o.c 

-0.600 
0.0 

o.c 
o.c 

0.0 
-0.102 
0.0 

0.599 
O.C 

o.c 

0.0 

0.469 

O.C 

0.0 
C.580 
-0.569 

0.0 

o.c 

0.0 

o.c 
o.c 

-0o?.46 
-Co  249 

O.C 

OoU 

0.0 

0.597 

0.0 

0.0 

o.c 
o.c 
o.c 

0.0 

o.c 
o.c 
o.c 


RATIO 

K 

L 

0.600 

10 

2 

C.  600 

11 

2 

0.600 

12 

2 

0.600 

13 

2 

0.600 

14 

2 

0.600 

9 

3 

0.60C 

10 

3 

0.600 

11 

3 

G.6C0 

12 

3 

0.600 

13 

3 

C.60C 

14 

3 

0.60C 

15 

3 

0.600 

7 

4 

0.600 

8 

4 

0.600 

9 

4 

0.600 

10 

4 

C.600 

11 

4 

0.600 

12 

4 

G.  6C0 

13 

4 

0.60C 

14 

4 

0.600 

5 

5 

0.600 

6 

5 

0.60  c 

7 

5 

Co  600 

8 

5 

0.600 

9 

c 

G.60L 

10 

5 

Co  6viO 

11 

5 

0.600 

12 

5 

0.600 

13 

5 

C.  60c 

14 

5 

0.600 

15 

5 

U.  600 

8 

6 

0.600 

9 

6 

0  e  60  0 

1C 

6 

.*:   AM;  . 

.  O  ■-  -    w 

1  i 

6 

0.6c0 

12 

6 

C.600 

13 

6 

0.  600 

14 

6 
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FINAL  VALUES  FOR  GRID  POINTS 


TEMP 

K 

I 

9.69 

1 

9.89 

2 

9.69 

3 

9.89 

4 

9.89 

5 

9.89 

6 

9.89 

7 

9.89 

8 

9.89 

9 

9.  64 

10 

8.79 

11 

11.05 

12 

12.88 

13 

11.46 

14 

9.89 

15 

9.89 

16 

9.89 

1 

2 

10.62 

2 

2 

10.62 

3 

2 

9.89 

4 

2 

9.89 

5 

2 

9.89 

6 

2 

9.89 

7 

2 

9.89 

6 

2 

12.45 

9 

2 

12.C4 

10 

2 

10.46 

11 

2 

ID. 66 

12 

2 

10.81 

13 

2 

8.96 

14 

2 

7.98 

i  c; 

2 

9.69 

1  6 

2 

9.89 

1 

3 

9.89 

2 

3 

9.89 

3 

3 

11.07 

4 

3 

12.06 

5 

3 

9.89 

6 

3 

1C.61 

7 

3 

12.22 

8 

3 

12.30 

9 

3 

10.89 

10 

3 

10.41 

11 

3 

10.34 

12 

3 

10.47 

13 

3 

10.39 

14 

3 

10.92 

15 

3 

12.90 

16 

3 

9.89 

1 

4 

9.89 

2 

4 

9.89 

3 

4 

9.89 

4 

4 

12.27 

5 

4 

11.43 

6 

4 

10.75 

7 

4 

11.15 

8 

4 

11.02 

9 

4 

9.06 

10 

4 

9.64 

11 

4 

9.92 

12 

4 

10.34 

13 

4 

9.84 

14 

4 

12.35 

15 

4 

9.89 

16 

4 
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9.89 

9.89 

9.89 

10.19 

7.51 

10.27 

11.43 

10.54 

10.26 

9.35 

9.55 

1C.38 

10.54 

10.45 

9.30 

6.45 

9.89 

9.89 

9.89 

10.01 

8.12 

10.72 

8.92 

10.39 

1C.10 

8.57 

11.01 

10.13 

11.02 

12.21 

6.03 

9.89 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 


5 

16 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 


5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
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MESH 
2 


MINIMUM 
6.000 


MAXIMUM 
14.000 


PROGRAM  OUTPUT 

MEAN  LAT 
37.000 


28  April   -   4  May 


ROTATION    PARAM    GRID 
47.00  1  1 


INT 
10 


INITIAL 
8. COO 


INCREMENT 
0.500 


POSITION    OF    THE    GRID    REFERENCE    POINT 
LATITUDE  LONGITUDE 

37  43c 200      124  6. COO 

OCEANOGRAPHIC    STATIONS 


EPSILON 
0.C01 


LATITUDE 

36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
36 
-  36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

36 

37 

37 

37 

37 

37 

37 

37 

36 

36 

36 

36 

37 

37 

37 


40c 
42c 
43o 
45. 
46. 
47  c 
49o 
51. 
52. 
54. 
55. 
56o 
54o 
53o 
51  o 
47c 
45c 
44c 
42c 
41c 
39 


50 

20 

80 

QO 

50 

5^ 

81 

40 

90 

10 

70 

20 

80 

50 

80 

10 

3n 

30 

,  6v 

20 

,40 


38.-20 


00 
60 
30 
30 
70 
70 
00 


36 

33 

31 

35 

39c 

42 

47 

50.00 

51.60 

53. 

54. 

55. 

55. 

0. 

3. 

4. 

2c 
31. 
33. 

Oc 

59. 
56. 
53. 
58. 
32. 
10 


3: 

30 

7  0 

90 

4C 

80 

20 

70 

40 

10 

10 

40 

30 

20 

40 

60 

CO 


10.00 


LONGITUDE 

121 

121 

121 

121 

121 

121 

121 

121 

121 

121 

121 

121 

121 

121 

121 

i2i 

121 

121 

121 

121 

121 

121 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 


DEPTH 


TEMP   SOUND  VELOC 


53. 7C 

54,  70 

54.  60 

54.30 

54.30 

54.30 

54.  5u 

54.30 

54.10 

54.40 

54. 

5°, 

59. 

59. 

59. 


58. 
^9c 
59. 
59. 
59. 


40 

30 
10 
40 
10 

40 

7  0 
10 

40 

20 
?r\ 

10 

60 
60 
00 

70 
20 
70 
70 
2C 


1, 

4, 

7. 

10, 

15. 

18, 

14. 

10.  _ 

9.70 

9.  20 

7.  50 

6.30 

a,  80 

14.30 

17.90 

19.40 

20.50 

21.80 

23.30 

23. 3C 

24.00 

28.60 

33.40 

40.50 

48.40 

45.80 
41.80 


10.00 


00 


10.00 
11.00 

11.00 

10.00 

9.00 
11.00 
9-  00 

11 

9.CC 

9.00 

9.00 

10  .CO 

10.00 

10.00 

1  o.oo 

10.00 

11.00 
10.0c 

11.00 

12.00 
10  c  CO 
11.00 
8.  00 
10.00 
8.00 
10.00 
10.00 
1 0  o  Q  0 
7.00 
8  •  C  0 
9.00 
11.00 
11.00 
10.00 
10.00 
8.00 
7.00 
8.00 
9.00 
9.00 
9.00 
11.00 
8.00 
9.00 
9.00 
10.00 
10.00 
9. CO 


10. 

10  o 

9. 
10. 
10. 

10 
10. 

11. 
11. 

He 

11. 

10. 
10. 

9. 
9n 


61 
51 
.76 
.28 
.27 
.03 


58 

03 

18 
rjg 

01 
63 

42 
85 
78 


60 
87 
08 
41 
63 
51 
47 
07 
\c 


9 

9. 
10 
10. 
!<*>. 
10 
10 

11 

10.  V 

10.38 

10.45 

9,98 
o 

9. 

9 
10 

9. 

9. 70 
10.3° 

9.33 

9. 

9. 

9. 

9. 

9. 

9. 

9. 
10 
10 


4P 
77 
14 
AC 


47 

39 

36 
20 

60 
60 

84 

2k 
46 


20 
9C 


10 

20 

90 

6C 

.80 


11.57 
12.45 
12.23 

12.1 


1490.00 

1469.40 

1487.00 

1488.80 

1488.60 

1487.80 

1489.80 

1491.30 

1491.90 

1491.50 

1491, 

1489, 

1489.20 

1487.30 

14  86.90 

1  A.  U  A  -  4 !_' 

I48  6. 30 

1^87.20 

14  8  3f 

1489, 

1489. 

1488, 

I486, 

1491.40 

1487.6 

148  8.50 

1488.90 

1487c 40 

1486.90 

1485.80 

1436.70 

1483.30 

1485.40 

1487.10 

148  8.70 

1486.50 

1486.70 

148  6.20 

1486.40 

1486.00 

1487.10 

1487.10 

1488.20 

1489.60 

1489.83 

1493.20 

I405. 80 

1475.20 

1494.70 


0.0 
0.0 

0.0 
0.0 

0.0 
0.0 
0.0 
0.0 

CO 

c .  0 

0.0 

0.0 

0.0 
0.0 

0.0 
0.0 

':,  e  KJ 
0  o  0 

0.0 

0.0 
c .  0 

0.0 
0.0 

0.0 
0.0 

0.0 

C  O    J 

0.0 
0.0 

0.0 
0.0 
CO 
0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

0.0 

0.0 
0  .  0 
O.o 
0  o  0 
0.0 
Co  0 
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122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

122 

123 

123 

123 

123 

122 

122 

122 

122 

122 

122 

122 


38. 5C 
35o20 
33.20 
31.30 
29.10 
26.70 
25.20 
27.50 
2^.60 
31.70 
33.  AC 
35.20 
37.  10 
38.70 
41.30 
43.10 
46.90 
49.90 
54.  50 
58.00 
3.30 
9.30 
15.40 
7.70 
59.70 
52.  50 
49.20 
48.90 
46.90 
44.  5C 
42.  10 


9.00 
10.00 
11.00 
10.00 
9.00 
10.00 
10.  oc 
9.00 
9.  00 
10.00 
8.00 
9.00 
9.00 
11.00 
11.00 
o.OO 
10.00 
9.00 
9.00 
9.00 
9.00 
8.00 
10.00 
11.00 
10.00 
10.00 
8. CO 
10.00 
9.00 
9.  r>0 
10.00 


11. 

18 

10. 

27 

10. 

57 

9. 

1C 

10. 

09 

9. 

49 

9. 

32 

10. 

05 

9. 

87 

9. 

51 

8. 

99 

9. 

32 

9. 

78 

9, 

52 

9« 

73 

10c 

59 

9e 

,75 

10, 

,46 

10, 

,91 

11. 

,04 

12.2 

11, 

,98 

10 

,61 

9 

,36 

9 

a65 

9 

.95 

9 

.92 

9 

.44 

8 

.93 

10 

.22 

9 

.82 

1491.80 

1489.30 

1490. 50 

1485.80 

1488.90 

1486.60 

1486.20 

1488.80 

1488.10 

1486.70 

1485.60 

1465.90 

1487.10 

1^86,80 

1487.40 

1490.60 

1487.40 

1489.70 

1491.30 

1491.40 

1495.30 

1494.30 

1490.3C 

148  5.40 

1^87.30 

1488.40 

1498.4" 

1  43  6.  4r> 

1494.60 

14  98.80 

1487.90 


CO 
0.0 
CO 
0.0 
0.0 
CO 
O.D 
0.0 
0.0 
0.0 

o.o 
o.o 

CO 
0.0 

o.o 

0.0 
CO 
CO 
0.0 
0.0 
0.0 
0.0 
OcO 
0.0 
0.0 
0.0 
0.0 
0.0 

U  e  0 
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FINAL  VALUES  FOR  GRID  POINTS 


TEMP 

K 

I 

9o70 

1 

9o70 

2 

9o70 

3 

9o70 

4 

9.70 

5 

9,70 

6 

9o70 

7 

9.70 

8 

13.42 

9 

lie  84 

10 

1^.32 

11 

8.79 

12 

7.62 

13 

8o09 

14 

8028 

15 

9o70 

16 

9o70 

1 

2 

9o70 

2 

2 

9o70 

3 

2 

9.70 

4 

2 

9.70 

5 

2 

12.77 

6 

2 

9.70 

7 

2 

9.7C 

8 

2 

12.83 

9 

2 

11.57 

10 

2 

10.46 

11 

2 

9.65 

12 

2 

9.32 

13 

2 

9,23 

14 

2 

11.75 

15 

n 

9  .  70 

16 

2 

9.70 

1 

3 

9.70 

2 

3 

6.10 

3 

3 

10.61 

4 

3 

11. c8 

5 

3 

10.78 

6 

3 

9.03 

7 

3 

12.23 

8 

3 

12.10 

9 

3 

11. C9 

10 

3 

10.28 

11 

3 

9.98 

12 

3 

10.  40 

13 

3 

10.38 

14 

3 

10.55 

15 

^ 

9.  51 

16 

3 

9.70 

1 

4 

9.70 

2 

4 

9.70 

3 

4 

9.36 

4 

4 

12.25 

5 

4 

9.57 

6 

4 

9.80 

7 

4 

10.29 

8 

4 

10.05 

9 

4 

10.42 

10 

4 

9.59 

11 

4 

9.  59 

12 

4 

9.81 

13 

4 

9.61 

14 

4 

10.54 

15 

4 

6.98 

16 

4 
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9.70 
9o70 
11.61 
9o65 
7.67 
9038 
10.20 
9.86 
9.92 
9.32 
9.47 
10.  14 
9.89 
9.86 
7.2C 
8.26 
9.70 
9.70 
9.73 
9.90 
8.26 
10.19 
9.76 
9.14 
10.46 
7.49 
10  .  63 
7.64 
11.01 
12.64 
8.10 
9.70 


1 
2 

3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 


5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

6 

6 

6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
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PROGRAM    FOP    OBJECTIVE    ANALYSIS    OF    CONSERVATIVE    PARAMETERS 
OF    THE    OCEAN 

REAL*8    TITLE 
fe0MMON/O6S/XO(120),Y0(120),OBSl(12C) 

2g8KJ8i)S8iU8 5x1  J8J!*?I  !  f| kj    10°' 

DIMENSION    CL(IO) ,LTG(3) ,TITLE(12) 
C    PK»R.8EglMSj5|SJm   °£ D|N FSorMfisiJ!L8ATBA?VTAR?RAB2E 

C    1    Srillfo-LlliV^i    MlcWISS^Ite"?5    THE 

c  RCTABffi  fl:SlJcSHfcSHJ.;teiftlW«in 

C  RICK  IS  THE MEAN  LATITUDE  FOR  THE  GRID 

WRITE( 6' III]    MESH, RMIN,RMAX, RICK 

giSBill^ll  AN^IPA^MI, INT, SINT, DELTA, EPSI 

ySfTil6,f430lANG,IPARfIM,INT,SINT.0ELTA,EPSI 
READ(5,292)  IL  AT  ,RMINIL,  ILONG,RMING 
WRITE  {  a  1  29C  ) 

WRITE < 6  I  293 )IL AT , RMINL, I  LONG, RM ING 
READ(5,280)XSC,YSC,NMR,MMR,NR 

THET=0o01745*ANG 
RICK=0o01745*RICK 

JSC  =  1 

READ*5,272)   (LATC I > ,0M INL < I > , LONG( I) , OMING(I> ,DEPTH( I 

2,OBl(I)70B2(I),OB3(I),I=ltNR) 

WRITE(6,270) 
WRITE  (6,271) 


2 

C  TO  0 


BrITEC*?*")  (  LAT< I ) ,OMINL( H , LONG! I > ,OMING< I > ,DEPTH< I 
fOBl(IJ»OB2U);CB3(I),I  =  lfNR)     RFFcRENrE  POINT  WITH 

ewafflWclM 

C    ROTATED,ANG,DtGREES    AND    TENTHS    CFPLuKto^T  NUM3ER    0F 

CC    S^ER^IO^       0      E    CH    S?£toS*88    ttlg    VERTICAL    GRID.FOR 
C    HORIZONTAL    GRID    MI     MUST    BE    ONEo 


WRITE     (6,321) 
WRITE(6,322) 
nn    iiQ     1=1*  NR 

403  IF(ILAT-LAT(J)  )404,406,^05 

404  LAT( J )=LAT ( J )-l  ^    _ 
CMINL( J)=0MINL(J)+60o0 

GO    TO    40  3 
4C5    LAT( J)=LAT( J )+l 

CMINL(J)=0MINL(J)-60o0 

406    IF( I  LONG-LONG ( J )  ) 407, 409, 40 8 
4C8    LCNG( J)=L0NG(J)+1     ^     . 
CMING( J)=OMING(J )-60o0 

GO    TO    4C6 
4C7    LONG( J) =LCNG( J )-!       „    „ 
CMING(J)=0MING(J)+60oC 

GO    TO    406 
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409  XO( J)=(RMING-OMING< J) )*C0S(RICK) 

Y0REL«X0(J)*SIN(tHET!+Y0(J)*C0S(THET) 
IF(MIoEQ.l)  GO  TO  411 
r  THF  VFRTICAL  GRID  IS  OBTAINED 

XO<J)=<(XORfU*2  +  YOPEL**2)**0. 51*1852.4 

Y0( J)=DEPTH(J) 
GO  TO  41«"< 

411  X0(J)=XGREL 
Y0(J)=YCREL 

410  IF(IPAR-2)412,414,416 

412  OBSK  J)=0B1(J) 
GO  TO  418 

414  OBSK  J)=0B2(  J) 
GO  TO  413 

418  WRITE(6;323)X0(J),Y0(J)t0BSKJ) 

419  CONTINUE 

2    XC=XSC/2oO 

YC=YSC/2oO  ,     ^    e 

XR=(XC**2+YC**2 )**0.5 

V)RITE(6,?49) 

WRITE(6,322) 

DO    15    L=1,MMR 

YG=YSC*L 
4    DO    15    K=1,NMR 

XG=XSC*K 

N  =  l 

DO    8    1=1, NR 

XOB=ABS(XO( I )-XG) 

IF(XOBoGT.XC)     GO    TO    8 

YOB=ABS (Y0( I )-YG) 

IF(YOBoGToYC)     GO    TO    8 

Z(N)=0BS1( I ) 

X  (  N )  =  X  0  (  i  J  -  X  G 

Y(N)=YO(  D-YG 
5    N=N  +  1 
8    CONTINUE 

IF(N-2)  1C, 12,14 
10    GVALHK,U=GoO 

12    AP=(X(1)**2+Y( l)**2)**0o5 

IF(AP.GTo7cC)     GO    TO    10 

GVAL1(K,L)  =  Z(1  ) 
IF(GVALl(K,L)aLToRMAX) 

GVALl(K,L)=CoO 

GO    TO    1 5 

14  N=N-1 
CALL    INTERP(N,XR,POLY) 
GVAL1 (K,L)=PCLY 
IF(GVALKK,L)oLToRMAX) 

15  WRITE(6!350()XG,YG,GV&L1(K,L) 

N1=NMR-1 


GO  TO  15 


GO  TO  15 


N2=NMR-2 

M1=MMR-1 

M2=MMR-2 
C  INITIAL  VALUE  FOR  U/K 
C  P  IS  THE  RATIO  OF  THE 

RAT 10=0 o 6 

WRITE  (6,319) 

DO  150  M=1,NJ1 

DO  150  L=1,MMR 

P=0o27 

RA  X=n  o  1 

RAY=C*( RATIO**2-RAX**2 )**0o5 . 
41  DO  15C  K=1,NMR  ,  , 

IF(GVALl(K,L).GToOol)  GO  TO  14' 

IF(MUEQol)  GO  TO  50 

IF(GVAL1(K,L+1 ) .LToOoCl)  GO 

IF(GVALl(K,L+2)  oLToOo?!).  GO 


TWO  COMPONENTS  OF  THE  VELOCITY  VECTO 


TO 
TO 


11C4 

11C4 
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GVALl(K,L)  =  PAY*YSC^(GVALl(KtL  +  l)-GVALKKTL+2))+2oO^ 

2IF(ABS(GVAL1(K,L))oLTbRMIN 
IF(ABS(GVALKK,L))oGT.RMAX) 


GO 
GO 


TO 
TO 


1102 
1 10  2 


1102 
1104 


GO  TO  149 
GVALKK,L)=OoO 


TO 

TO 


149 
1106 


IF(GVALl(K,L-ljoLToOo01)GO 
IF(GVALKK,L-2)eLT,C»01)GO 

g5ALlCK,L)-RAY*YSC*(GVALlCK,L-2l-GVALl(KfL-in+2.0* 


2GVALl(KfL-l)-GVALl(K,L-2) 
1F(ABS(GVAL1  KtL  -LT.RMI 
IF(ABS(GVAL1(K,L)  UoT.RMA 


v|  )  GO 

AX)  GO 


TO 
TO 


110  5 
1105 


GO    TO    149 
1105    GVALl(K,L)=OoO 
1106 


56  ^mkHfiVSWcW^t-H?! y^A^+»  l/2-°-R*y*vsc* 

2(GVAL] (KiL-1 )-GVALi(K,L+l)     /4of 


1108 
50 
51 
54 
42 
80 
83 
81 


C-,0    TD    149 

GVALK  K,L)=OoO 

GU   TO    149 

N=0 

IF(NMR-(K+1))    51*54 

IfiLoEQ-riHRJ    GO   Tf) 


,42 
18  5 


TO 


GO 
1  f  .  L 
IF(L 
IF(L 
£F(L 


l£Q«MMR)    GO 
EQ.li    GO     «0 
GT.M1)    GO    TO     lb 
EQoi)    GO    TO    110 


0    01 
11C 


IF(GVAL1(K+1  iL+D.LToO 
IF(GVALl(K+2»L)oLTal;ol 
I  F  (  G  V  A  L  1  (  K  + 1  » L  - 1  )  o  L  !  o  u 


,1IG0    TO    81 
)    GO    TO    81 
ol)G0    TO    81 


N  =  7 


GO 
GU 


TO 
TO 


1.0r8 
IOCS 


1C08 
ICO 


103 
104 


IF(GVALl(K+l,LULTo^oi)GO    TO    103 
IF(NoLT.7)     GO    TO    IOC 

CALL    MOLE(MMfK,L,XSC,YSC,OIV) 

2![f1ab^(gvali(k1l|j.gt.rmaxj 

IF ( A3S ( GVAL1 ( K t L ) ) o  LT. RMI N) 

GO    TO    149 
GVALl(K,L)=CoO 

IF(GVALKKtL-l  )•  LTo 
IF(LoLQoMMR)    GO    TO 
Ir(GVALl(K,L+l JoLTo^ol) 

N=7 

GG    TO    104 

riF?K0LT„2)    GO,  TO    132 
IriGVALl  (K-liL).LI  o'^.l) 
IF     (N.LT.7)    GO    TO    108 

MM- 9 

c 


>ol) 

104 


GO    TO 
GO    TO 


103 
103 


GO    TO    110 


CALL    MOLE(MM,K,L,XSC,YSC,DIV) 
DIV=(GVAL1(K+1,L)-GVAL1(K-1,L)) 
K , L+l ) -GV AL1 ( K , L- 1 )    *RAY/ ( 2 o 0-' 
GVAL1(K,L)=GVAL1<K,L|-0IV 
If { ABS (GVAL1 ( K  »L    ) oGT0 RMAX) 
I F ; ABS ( GVAL1 ( K »  L ) \  . L • • RrfIN ) 


)*RAX/(2oO*XSC)+(GVALl( 
<YSC) 


GO 
GO 


TO 
TO 


1009 

1009 


10  09 
1^8 


llr> 
111 


-0.1 
GO 


Tn 


GO    TQ    149 

G^ALH  K,LJ 

IF [L.Li .2J 

»F  'GVALl (K-lt L-l 

GO    ro   ill 
N-0  .    _., 

Ir.KoLToJ*     GO     rO 
,     v GVALl (K-2t L) o 


ilO 

oLT, 


Ool)    GO    TO    110 


LT 


3  2 


1) 


TO    114 


75 


112 


I:-  :Lo  EQ<  W1R) 
ir \ GV«Li  i K 


GO 
ltL+1 


m  li* 

TO  113 


114 


1004 

114 

113 
115 


132 


GO  TO  132 


GO  TO  134 


ICO  3 
134 

136 


1)  GO  TO 
IF(N«LT«7) 

TALL  MOLE(MM,K,L,XSC,YSC,DIV) 

8fe!Ell^tl«:tl:M«?K*LiiK-i.Lni-«»LHK-i. 

GO  TO  149 
GVAL1(K,L)=0«0 

GO  TO   113 

N=0 

GO  TO  115 

IF(LoEQoMl)  GO  TO  75 
IF  N.LT.7)  GO  TO  132 
IF(GVAL1(K-2,L+2).L i.o.l) 

N  =  N+8 

IF1L.GT.M2)  GO  TO  136 
IF(GVALl(K,L+2 J.LT.Q.l) 
IFCN.LT.15)  GO  TO  136 

TALL  MOLE  WtK,L,XSC,YSC,DIV) 
SvttlJKTL  =G§Atl  KfLI-RATIO»D|V 

GVALKK, L)=bVALl (K,L) 

Sf(ABS(GVAL1(K,L)JoGT.RMAX) 

IF(ABS(GVALl(KtU)oLToRMIN) 

GO    TO    149 

GVALl(K,L)=OoO 

IF(K.EQ.NMR)    GO    TO    75 

GO    TO    149 

N=0 

IF(L.GT.Ml)    GO    TO    75 

IF(K.GT.Nl)     GO    TO    75 

IF(GVAL1(K+1,L+1 loLTo' 

IF(LoGT.M2)    GO    TO    1*8 

IF(K.GT.N2)^GO    TO    148 

IF(GVAL1(K  +  2,U.LT.«:.1)G0    TO 

IF(GVALl(K+2,L+2).LT.i'.l)     GO 

MM=1 
CALL    MO 
GVALKK 
GVALKK. 
L))/(2.0*XR) 


-RATIO*  (  GVALKK,  L  +  2)-GVALKK-2, 


GO 
GO 


TO 
TO 


1003 
IOC  3 


DGO    TO    75 


148 
TO 


148 


II  E(MM,K,L,XSC,YSC,DIV) 
tL    =GVAL.l(K,L)'RATIO*DIV 

,L)=GVAL11K,L)    -RATIO*(G 


VALKK,L  +  2>-GVALl(K+2, 


GO 
GO 


TO 
TO 


10^1 

ioAi 


1001 
148 


I  F ( AB  S ( GVAL1 ( K , L ) )  « GT. RMAX) 
IF(ABS(GVALl(K,L))oLT.RMIN) 

GO    TO    149 
GVALKK,  L)=C.O 
IF(KoGToNl)     GO    TO    75 
IF(L.GT.Kl)    GO    TO    75 
IF(KcLTo2)    GO    TO    75^ 
IF(GVAL1(K,L  +  1  )oLT.'.. 
1F(GVAL1(K-1,L+1).LT, 

CALL    MOLE(MV,K,L,XSC,YSC,DIV) 

8teB:h!3fel:l«:tl-sa5S«»ALicK-i.L*i.-*»*Li«Kii. 

L+1))/(2.0*XSC) 


1) 


GO    TO    75 
1)     GO    TO    75 


GO 
GO 


TO 
TO 


1002 
75 

82 


185 


IF(ABS(  GVALKK,  LoGToRMAX) 
IF(ABS(GVALHK,L))oLT.RMIN» 

GO    TO    149 

GVALKK, L)=0o0 

IF(LoLT.3)     GO    T0    149 

IFCK.GT.N2)     GO    TO    185 

IF(GVALKK  +  2,L  )oLToOol)     GO    TO     185 

IF(KoEO.NMR)    GO    TO    85 
IF(GVAL1(K+1,L-1  )oLT.0.1)GO 

N=N+2 


1002 
1002 


TO    85 


76 


85 


LT.0.1)  GO  TO  149 


GO  TO  91 


1C07 
91 


94 


IF(GVALl(K,L-2) 

TfTkUeoHZ)    GO  TO  91 
IF(GVAL1(K+2»L-2).LT.0.H 

IF(N.LT. 7)     GO    TO    91 

GO    TO    149 

^HftflrKVtf.O.ll    GO    TO    1« 

IF(GVAL1  (K-lfL-1  ).L  t  ou  x  i 

?F(GVALl(K,t-l).LT.0.1)    GO    TO   96 

IF(NoLTol4)    GO    TO    96 

CALL    MOLEW,K,LjXSC,YSC,DIV) 
L-l).)/(2-0*XSC) 


GO 
GO 


TO 
TO 


10^6 

10°  6 


1006 
96 


1005 
140 


IF(ABS(GVALl(K,L))cLToRMIN) 

GO    TO    149 

GVALl(K,L)=(,oO  n    ,.     G0    T0    i^g 

rVALl(K?Ly=GVALllK.ti-i       r\  h  i  -  <j       ^ 

2IFiAis(GVALl(KtL))oGT0^AX) 
1f(ABS(GVAL1(K,L)).LToRMIN) 

GO  TO  149 
GVALl(K,L)=CoO 

?F(m.EQ!l>GO  TO  142 

iSU  lift J -i StftSitT"*'"" ' 

2(YSC^(GVAL1(K,L-1)-^VAL1(K,L+D  I 
GO  TO  154 


GO 
GO 


TO 
TO 


TO 
TO 


100  5 


149 
149 


142  TF(GVALl(K,L-2|oLTo 


!ol)  GO  TO  149 
:,1)G0  TO  1-+9 


IF(GVAL1 (K-2»L ) . LT • 

2^lS?j^g^H^;?^SL,+GVALm-z'u+GViL1,K' 

2GVALKK?L-1)+GVAL1(K»L    2))/HKb 


rj05,\  ))    GO    TO    149 
;GEcRATI0)G3    TO    1021 


1021 


1C22 


10  23 

1C24 

153 


IF(ABS(ROX 

IF(ABS(ROX) 
RAX=ROX 
GO    TO    1022 
RAX=RATIO 
RAY=0.0 

b2(R°TI»S2-RAX**21**0.5 

IF(RAXoGT.OoO)GO    TO    1023 

RAY=-C*B 

GO    TO    1024 

IfTm^LT.KJI)    GO   TO    157 

IF(ABS(RAY)oLTo|PSI )  GO 
IF(A3S(RAX)oLToEPSl)  GO 
p=RAY/RAX 


TO 
TO 


1053 

10  54 


77 


10  5: 


1C55 
154 
149 
157 

150 


LToEPSI )    GO    TO    1055 


GO    TO    154    , 
IF(ABS<  RAX)„ 
P=EPSI/RAX 
GO    TO    154 
!C54    p=RAY/EPSI 
GO    TO    154 

WRITE<6,323)     RAX,RAY,RATI 0, K , L 
IF(MoLToMJl)    GO    TO    157 
IF(GVAL1(K,LUGT»RMIN) 

GVALl(KtL)=OoC 
CONTINUE 

JSCoEQoMESH) 


GO    TO    150 


GO    TO    170 


MMR 
NMR 


0.1)    GO    TO    162 


IF( 

N=l 

SUM=0oC 

DO    162    L=lt 

YG=YSC*L 

00    162    K=! 

XG=XSC-K 

IF(GVALKKtL)oLT. 

03S1(NP+N)=GVAL1 (K,L) 

Y0(NR+N)=YG 

XC(NR+N)=XG 

N=N+1 

SUM=SUM+GVAL1( K,L) 
162  CONTINUE       T  ,  ki  ,  .  . 

AVER=SUM/(FL0AT(N-1)) 

NR=NR+N 

JSC=JSC+1 

NMR=2*NMR 

MMR=2*MMR 

XSC=XSC/2.C 

YSC=YSC/2,0 

MJ1=3*MJ1 

WRITE (6, 280 

GO  I U  ? 
170  DO  173  L  =  l, 

I3=MMR+1 

I=I3-L 

YG=YSC";L 

DO  173  K=1,I\MR 

IF(GVAL1(K,L).GT.0.1) 

GVAL1 ( K, L)=AVER 

172  XGA=XG*r0S(THET)+YG*SIN(THET) , 
1    YGA--XG*SIN(THET)+YG*COS(THET) 

XGA=RMING-XGA 
YGA=RMINL+YGA 
G(Kt I)=GVAL1(K,L ) 

173  CONTINUE 
WRITE(6,255) 
WRITE(6,256) 
WRITE(6,259) 
IDR=INT-1 
DO  171  1=1, IDR 
CL( 1)=SINT 

171  CL( I+1)=CL( I )+DELTA 
LTG(l)=IS.ECo2 
LTG(2)=JS.E0o2 
LTG(3)=KSoEQo2 


,  j  )    X  SC .  Y  S C  ,  NM  R ,  MMR  , N R 

MMR 


GO  TO  172 


((GVALl(K,L),K,L,K=l,NMR)tL=l,MMR) 


204 
255 
256 
258 
2  59 
27fi 

271 

272 
273 
274 


*CCNTUR(GVAL1,NMR,MMR,5C,CL,INT,TITLE,9,21,LTG) 


24X^    [NAL  VALUES  FOR  GRID  POINTS,/) 
3r  X , 4HT  E MP , 5  X , 1 HK , 5  X , 1 HL ) 


CALL 
F0RMAT(F10o4,F15«4,fl0.2) 

FORMAT     CO'     "" 
FORMAT ( 'O1 , 
F0RMAT(2X,        F15o2,3I5 

PoRMaV(?0«;40X?22HOCEANOGRAPHIC    STATIONS,/)     _mf=OTH    5X 
F0RMAT(«'";,;2^X|8HLA-TITUDE,7X,9HL0NGITUDE,8X,5HDEPTH,5X 

?    4HTFMP,2X,11HSCUND    VELOC,/) 
"  FORMAT (  I5,F10« 2,I5,5F10o2) 

FQRMAT(12X,3F15«4> 

FORMAT (20X, I  5 ,F10. 2 , 15 , 5F10.2 ) 
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280 
284 

2  90 
291 
292 
293 
319 
321 
322 
323 
324 
325 
326 
349 

3  50 
401 
40 


FORMAT!  2FKo0t3I5) 

^RMATMC«i20X:36HP0SITI0N  OF  THE  GRID  REFERENCE  POINT 

FORMAT( 26X: 8HL ATI TUDE , 5X ,9HL0NG I TUDE , / ) 

FO°MAT ( I5,F10«0,  15 t  P1^«L ) 

FORMAT!  20X;iHX,15X,lHYaOXt9HPARAMETER,/) 
FORMAT! 10X,3F15« 3, 215) 

FORMAT!  if5^2rX°,4HMESHf2X,7HMINIMUMf3X,7HMAXIMUM,2X,8HM 

FORMAtI'i'UFx^^I^ERPOLATED    VALUES    FOP    THE    GRID,/) 

FORMATtl2X,3F15.5) 
F0RMAT!F10o2,3I5,3F10.3) 


4  30 

888 
890 

898  F0RMAT!6A8) 

899  F0RMAT!3I5) 
STOP 

END 


•:x.o 


SUBROUTINE  INT ERP ( N , XP , POLY ) 
PROGRAM  FOR  LEAST  SQUARE  SURFACE 

COMMON/ PAR AM/X I 2U) , Y! 20)  ,  Z I  20 ) 

DIMENSION    C(6,6) ,B!6) 

IF(NoLTc3)     GO    TO    2^ 

M=6 

DO    2    1=1,6 

R(T )=0.0 

DO    2    K=l,6 
2    C!K,I)=OoO 
C(ltl)=N 

DO    5    J=ltN       ,     Kti  1V 
C(1,2)=C(1,2    +X  U 
C!1,3)=C!1,3    +Y    J 
C!lr4)=C!l,4)+X    J) 
CU,5)=CU,5)+XU   *Y!J) 
C!1,6)=C!1,6)+Y!J**2 
C(2,4)=C!2,4)+X! J)**3 
Clil?   =C!2,5)+Y!J)*X!J)**2 
C    2    6    =  C!2    6)+X(J)*Y!J)**2 
C(3,6)=C!3,6)+Y!J    **3 
C!4,4)=C!4t4)+X!  J)**4 
C    4,5 =  C    4 ,5)+Y!J)*X! J)**3#     v  .^ 
C    4    6    =  C    J    6+1X1 J)**  2)*(YU)**2) 
C    5   6   =C!5;6)+X!J)*Y!J)**3 
C(6,6>=C!6,6)+Y!J)**4 

B!1)=B!1)+Z!J) 
B!2)=B!2)+Z!J)*X    J 
B(3)=B!3)+Z(J)-Y(J 
B!4)=B!4)+Z!J)*X   J)**2 
B(5)=B(5)+Z(J)-X(J)*Y    J) 
5    B(6)=B!6)+ZtJ)*Y!J)**2 
C!5,5)=C<4,6) 
C(3,5)=C(2,6) 
C(3,4)=C(2,5) 
C(3,3)=C!1,6) 
C(2,2)=C!1,4) 
C(2,3)=C<1,5) 
DO    7    1=1,6 


TO    1 


INTERPOLATION    UP    TO    SECON 


DO    7 

K= 

=  1,6 

C!K, 

I)  = 

=C!I  1 

,K) 

IF!N 

oGl 

.5) 

GO 

N  =  4 

DO    9 

1  = 

=1,6 

C!4, 

I)  = 

=  C(5, 

,  I) 
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9    C(  It4)=CUi5) 

10  CALL=SIM0  (C,B,M,KS) 
PCLY=B(1) 

-MyHll»fiBi|.S.».8»-'«»-"»»' 

2UX(l)-X(2ni(Y(l)-Y(2))) 

24  P0LY°(2Z(1I+  1(2)1/2.0 

||  W;}x  1)**2+Y(  1)**2I**0.5 
CS=1.5703*XY/XR 

POLY=Z<  D-CM 
42  RETURN 
END 

SUBROUTINE  HOLE  <MM,K,L , XSC.YSC.DIV) 

COHMON/GRID/y(50,50) 

FAC1=(XSC/YSC)**2 

,     U?J?,ii'^U^^;^Ii-Ui^,L"2)-U(K  +  2,t»    -U<K+2,L+2> 

1  HlV^(U(Kii,L;il-U(K->2,L>2U/XR 

2  §?K";°)i-FAC2*(U(K+l,L+l)-2.0*U<K,L+l)*U(K-l.L+ll)*2.C* 

,    R?KTPl-4    0*U(K-l,L+l)-U(K,L+2)-U(K-2,LI-U(K-2,L+2l 

3  ^<u7k-1*L«)HJ<K-2,L+2!>/XR 

4  §?kTS,15FaC1*(UIK-1,L+1»-2.0*U(K-1,L1*0(K-1.L-1I»«.C* 

^V^U.k^L^uiklz.D./XSC 

5  S?kTl)".0*U<K-1,L-1»-U(K-2.LI-U«K.L-21-U<K-2.L-2I 
DIV=(U(K-l,L-l)-U(K-2,L-2))/XR 

6S?kTC)"f»C2*IO(K*1.L-II-2.0*U(K.L-1)««K-1,L-1H+2.0* 

^V^fuT^rL-n-ulK'.L-ai./YSC 

7  8?KTe.it0*U«Ktl,iV-l.-U(Kt2,L»-U(K,L-2»-U.K+2,L-2, 

8  8?KTLli-FACl*CUCK+l,L+l)-2.0*UCK+l.L)+U(K*l,L-ll)+2.0* 

2oW5tftlfeilV£?-S!  Wz.Ll  1  /XSC 

9  S?K^)^U(K+llL)+FAClMU(K,L+l)+U(K,L-l))-U(K-ltL))/( 

2    2«,**FAC1  +  2«C) 
10    RETURN 
END 
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